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1. Introduction 



The lightest neutralino is one of the strongest candidates to explain dark matter. It natu- 
rally fits the definition of a weakly interacting massive particle (WIMP) and is the lightest 
supersymmetric particle in many models in the literature. Furthermore, the lightest super- 
symmetric particle (LSP) is stable in most viable models involving super symmetry. These 
properties allow the neutralino to provide a significant fraction of critical density to the 
universe. Throughout this paper, we assume that the lightest neutralino is the LSP and 
that it is the dark matter candidate whose relic density we will be calculating. 

Measurements of the dark matter relic density are becoming a powerful constraint on 
all extensions of the standard model of particle physics. Recent improvements in these 
measurements allow dark matter densities to be even more discriminating [Qj. In the most 
studied extension, mSUGRA, only a few viable regions of parameter space are left [g, 
Thus, it is becoming increasingly important to have accurate calculations of the dark matter 
density. 

In this paper, we present analytic cross sections for all tree level neutralino-neutralino 
coannihilations |||, [j| in the MSSM for general neutralinos. Since our calculations allow 
the possibility that the two interacting neutralinos are not identical, our results include 
and generalize previously published expressions. In an effort to standardize the notation 
of neutralino cross section calculations, and for ease of comparison, we closely follow the 
notation of ||. We also make no simplifying assumptions about a given MSSM model 
beyond neglecting most of the CP-violating phases in the SUSY sector. That is, we make no 
additional assumptions about the neutralino, chargino, or sfermion sectors. We also include 
all s-channel widths. Thus, these results provide a step forward towards a complete analytic 
calculation of the neutralino relic density without using the normal velocity expansion. 
Additionally, these results are valid regardless of proximity to resonances and final state 
thresholds. However, coannihilations between the neutralino and a chargino, a neutralino 
and a sfermion, and two charginos can also be important in relic density calculations. 
Results for squared matrix elements involving stau and stop coannihilation can be found 
elsewhere in the literature ^, |J . Complete cross section expressions involving charginos 
will be presented by the authors in the future. 

The calculated expressions are voluminous. The size of the expressions and the varying 
conventions in the literature forbid most direct analytical checks. However, as we have 
taken steps to follow most conventions of ||, we are able to check our answers in the limit 
of the two initial state neutralinos being identical. Our results match exactly. 

In addition to including the standard neutralino relic density contributions, these ex- 
pressions also include contributions which are important when the mass of the next-to- 
lightest neutralino is approximately degenerate with the lightest neutralino. In mSUGRA, 
this occurs in the higgsino-dominated region of parameter space characterized by large 
mo |1C]. As too is increased, the /i parameter eventually falls, leading to high higgsino 
content. If rriQ is increased even more, it quickly becomes impossible to achieve correct 
electroweak symmetry breaking. This difficulty in acquiring neutralino mass degeneracy 
in mSUGRA is directly related to having universal gaugino masses at the SUSY-breaking 
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scale. In contrast to mSUGRA, many models in the literature, including models derived 
from the heterotic string [11, 12, 12, 14] and the Type-I superstring [15], generically have 
non- universal gaugino masses. Non- universal gaugino masses allow for larger sections of 
parameter space where the two lightest neutralinos are almost degenerate. The dark matter 
prospects of non-universal gaugino masses have been studied in Q, where it was pointed 
out that such models compare quite favorably with mSUGRA in terms of both motivation 
and neutralino relic density. 

The outline of the paper is as follows. In Section 2 we review the relic density calcu- 
lation. Here, again, we follow the format of Q, which uses the simplifying work of [jO], 17] 
to motivate the format of the cross sections. In Section 3 we present the cross section ex- 
pressions for all tree level coannihilations of two neutralinos. These expressions are exact 
and include integration over the final state angular distribution. Since they are exact, they 
can be used even when the masses of the two neutralinos are far from degenerate. This 
allows the expressions to be used for purposes in addition to relic density calculations. We 
conclude in Section 4. 



2. Review of the Relic Density Calculation 



The calculation of a relic density from thermal considerations is standard [21]. One starts 
with a theory that contains a particle expected to have a significant relic density. One 
wishes to calculate the particle's contribution to the total matter density of the universe, 

= Pxl 1 Ptot = m x n xl Ptot- This number is then compared to the current experimental 
bounds on a particle's relic density. The limit for the dark matter density is generally taken 
to be 0.1 < Q x h 2 < 0.3 though more aggressive limits have been promoted M. 

The initial number density of the neutralino is assumed to be determined because the 
particle is in thermal equilibrium. When the particle begins in thermal equilibrium with its 
surroundings, interactions that create neutralinos usually happen as frequently as reverse 
interactions which destroy neutralinos. Once the temperature drops below T ~ m x , most 
particles no longer have sufficient energy to create neutralinos. Now neutralinos can only 
annihilate, and these annihilations occur until about the time when the Hubble expansion 
parameter b ecomes larger than the annihilation rate, H ^ T ann . When expansion dwarfs 
annihilation, neutralinos are being separated apart from each other too quickly for annihi- 
lation to be efficient. This happens at the freeze-out temperature, usually at Tp ~ m x /20. 

In most relic density calculations, the only interaction cross sections that need to be 
calculated are annihilations of the type XX ~^ X where X is any final state involving only 
standard model particles. However, there are scenarios in which other particles in the 
thermal bath have important effects on the evolution of the neutralino relic density. Such 
a particle, called a coannihilator, usually must have direct interactions with the neutralino 
and also must be very nearly degenerate in mass. When this occurs, the neutralino and 
all relevant coannihilators form a coupled system. In this section, we will denote particles 
belonging to that coupled system by Xi- Now all interactions involving particles in this 
coupled system come into play, including XiXj ~^ X, XiX — ► XjY , and x% ~^ XjX. Here 
both X and Y denote states including standard model particles. Decays once again enter 
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the calculation because the coannihilators are generally not stable and eventually decay 
into the lightest neutralino. 

The effect of this coupled system on the relic density of neutralinos can best be under- 
stood through a consideration of the freeze-out temperature. Because these coannihilators 
have essentially the same mass as the neutralino and share common interactions, they re- 
main in thermal equilibrium for about as long as the neutralinos. The neutralino can thus 
interact through more channels with more particles than in the usual single-species anni- 
hilation scenario. This raises the effective annihilation rate, so the determining equation 
for freeze-out (H > T ann ) is satisfied at a later time and a lower temperature. The neu- 
tralino is nonrelativistic during freeze-out, so its equilibrium number density is decreasing 
exponentially with decreasing temperature. Since coannihilation keeps the neutralino in 
equilibrium down to a lower temperature, the number density at freeze-out is lower. This 
lowering of the number density is the most relevant effect of coannihilation Q and can 
have drastic consequences for the viable parameter space of many supersymmetric models. 
In the MSSM, many particles can serve as the coannihilator: the lightest stop, the lightest 
stau, any of the other neutralinos, or either of the charginos. The cases of the stop and 
stau have been studied extensively in the literature Jj], [|. Here we study the case of 
coannihilation between the lightest neutralino and another neutralino, reserving presen- 
tation of the results for coannihilation between the lightest neutralino and a chargino for 
future work. 

For the case without coannihilations, evolution of the particle number density happens 
in accordance with the single species Boltzmann equation: 



dn 
~dt 



-3Hn — (av) 



n 



(n 



(2.1) 



The number density is reduced by Hubble expansion and annihilations of the relic particle. 
The number density of the relic is also affected by its own decays and inverse annihilation 
processes. In most cases, including the versions of the MSSM studied here, the relic particle 
is assumed to be stable, so relic decay is neglected. In the above expression, we have also 
assumed T invariance to relate annihilation and inverse annihilation processes. As noted 
earlier, the thermally averaged cross section (av) involves only processes such as XiXi ~^ X. 

In the presence of coannihilators, the Boltzmann equation becomes more complicated 
(Equation (27) from |] 



drii 
~dT 



N 



3Hm - ^(crijVij} (riiUj - n^nf 



Y [r^ (m - nf) - Tji (rij - nf 



(2.2) 



Here i and j can denote the relic particle or any possible coannihilators. The velocity Vij is 
the relative velocity between particle i and particle j. The first cross section aij is a sum 
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over final states of processes such as XiXj ~^ X. a' Xi j is a sum over standard model final 
state particles of XiX — > XjY- The V terms describe possible decays of the coannihilators 
Xi -> XjX. 

One uses the different stability properties of the relic particle and the coannihilators to 
simplify the expression. As stated above, the relic particle is generally assumed to be stable 
(usually from R-parity in the MSSM). The coannihilators' decay rates are also generally 
quite fast compared to the age of the universe. In the MSSM, R-parity also guarantees 
that each coannihilator will eventually decay into standard model particles and the lightest 
neutralino (being the LSP). Since the LSP is the relic particle that we are interested in, we 
can sum the Boltzmann equation over the index i to get an equation for the final number 
density of the relic particle GJ. Now we use n = YliLi n «- The terms in the second and 
third lines now cancel, so we have: 



dn N 

= - 3Hn - ^ (vijVij) {rurij - n^nyj (2.3) 

To a very good approximation, one can use the usual single species Boltzmann equa- 
tion for the case of coannihilations if the following replacement is made for the thermally 
averaged cross section: 



r> eq n eq 

The thermally averaged cross section deserves careful treatment. For scenarios involv- 
ing coannihilations, the expression for the thermally averaged cross section begins as a six 
dimensional integral, seven dimensional including the integration over the final state angle. 
This integral has been conveniently put into the form of a one-dimensional integral over 
the total squared center of mass energy p 



M = 7-TT— ^ , " "* 2 " ^ (2-5) 



8mjT 



Here (Tjj (s) is the sum of all cross sections for any process XiXj ~^ X and @f (s, mj, rnj) is 
a function related to the incoming and outgoing momenta: 



(3f (s,mi,mj) = ~\J{ s ~( m i + m j) 2 ) [ s ~ ( m i ~ m jY) ( 2 - 6 ) 

Consideration of 2-body final states is sufficient for almost all relic density calculations, 
that is what we consider in this paper. Part of the integrand can be written in the form: 



-fy (s, m^rrij) a {j (s) = — \c 6[s - (m fl + m /a ) 2 j f (s, m fl , m h ) (s) 

hh 

(2.7) 



- 4 - 



Here c is a color factor, equal to three if the final state is a quark-antiquark pair, equal to 
one otherwise. o(^s — (m^ + mf 2 ) 2 ^j is the usual step function to enforce the hard lower 
limit on the squared center of mass energy. The cj functions are the squared amplitudes 
averaged over the final state angle: 



*fih 0) = \f dcosOcM \A(ij - /i/ 2 )| 2 (2.8) 
Computer codes exist which numerically perform this final state angular integration. Some 



op- 



codes include only a subset of all possible coannihilation channels [18, 19] while the pro- 



gram micrOMEGAs includes all relevant coannihilation channels [20]. We have performed 
the integration analytically, thus reducing the number of necessary numerical integrations 
and allowing for easier analysis of relic density properties. The next section lists these 
u)ij-,f x f 2 (s) functions for all possible 2-body final states from neutralino (co)annihilation. 

Once the effective thermally averaged cross section has been determined as a function 
of temperature, it is used to iteratively calculate the freeze-out temperature. In practice, 
the dimensionless inverse freeze-out temperature xf = m x /Tp is calculated using: 

x F = In ( ^9m Pl m x {av) \ 

Here mpi is the Planck mass, g is the total number of degrees of freedom of the x particle 
(spin, color, etc.), g* is the total number of effective relativistic degrees of freedom at freeze- 
out, and the thermally averaged cross section is evaluated at the freeze-out temperature. 
For most cases in the MSSM, xf — 20. 

When the freeze-out temperature has been determined, the total (co) annihilation de- 
pletion of the neutralino number density can be calculated by integrating the thermally 
averaged cross section from freeze-out to the present temperature (essentially T = 0). The 
thermally averaged cross section appears in the formula for the relic density in the form of 
the integral J (xf) = f^i&v) x~ 2 dx as 



7T h 2 so 1 

o m Pi (g*s/gl /2 ) J(x F ) 



n xh 2 = 40 A /^ff 7 : (2.10) 



which is commonly given as 



n , 2 1.07 x lO 9 ^^- 1 

n x h 2 = — — — . (2.ii) 

gj m pi J(x F ) 

This is the expression that one compares with the relic density limit of 0.1 < Qh 2 < 0.3. 
Once a freeze out temperature has been determined, one must still integrate over x to 
determine the relic density. In practice, these integrations are usually done numerically. 
Previously, it was also necessary to integrate over the final state angular distribution. The 
results listed in the next section, along with those of H 0, have already performed this 
angular integration. Our results also represent a significant step towards a full tabulation 
of all analytic cross sections relevant in the accurate determination of a relic neutralino 
density. 
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3. Exact Cross Section Expressions 

In this section we list all possible Q functions for two neutralinos (co)annihilating. All of 
these processes involve exchange via s, t and u-channels. We use the following definitions 
for the sum and difference of the neutralino masses: 



a = mo + mo (3.1) 

Xi Xj 

5 = m v o — m v o 

Xi Xj 



We refer the reader to || for the definitions of the % and 3^ functions. Since we 
are studying a generalization of the processes studied in |J, their definitions for particle 
couplings are also used here. Furthermore, we suppress most of the functional dependence 
of the T functions. The suppressed dependencies are apparent from the % and ^ functions. 
We also refer the reader to || for an explicit listing. However, there are some places where 
the difference between the two neutralino masses comes into play. If the function F (listed 
in ||) is defined in terms of the incoming (p) and outgoing (k) CM momenta as F = 2 hp, 
then the only important difference is in the function D: 

n / \ Xi + x 2 yi + y 2 s {xi- x 2 ) {yi - y 2 ) ,„ oA 

D(s,x 1 ,x 2 ,y 1 ,y 2 ) = + — — (3.2) 

Here x\ and x 2 are the incoming squared particle masses and y\ and y 2 are the outgoing 
squared particle masses. 



Contributions to u> come from s-channel Higgs boson exchange, t- and u-channel neutralino 
exchange and cross terms 

~(h,H) ,~(x°) ,~(h,H- x °) ,„_v 
"xix^Hh = UJ Hh + u Hh +UJ hH ■ ( 3 - 3 ) 

S-channel CP-even Higgs bosons (h,H): 

{s - a 2 ) ; (3.4) 



r ,(h,H) _ i 

r=h,H 

T- and U-channel neutralino: 



(jrhH QXiXjl- - 



s — ml + im r T r 



~,(x°) 



Hh 



E 

fcj=i 



m x2 m x? J fc^ +m x° k J kl H 



(3.5) 



where 
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rhH 
1 kl 



jhH _ 
J kl — 



K, 



hH 
kl 



and 



s 2 -s(a 2 + 6 2 ) + a 2 5 2 



[C r Hh % - C y Hh y Q ] (s, a, 6, ml m 2 H , m 2 o, m 2 o) , (3.6) 



1 

4s 



-8C^ S ar 1 + GS {0) T 



- 4C Hh sa (yi + {s- a 2 ) y )) (s,a,6,m 2 h ,m 2 H ,m 2 x o,m 2 x o) , 



1 

16a 



LbC Hh Sl 2 + Lr Hh ll+Lr Hh JO 

+ lQCy h sy 2 - 8 (c% h - 2C y H i) o5 (m 2 H - m 2 ) y x 



(3.7) 



(3.8) 



G 



J,T(0) 
Hh 



8C^S{m 2 H -m 2 ) (s - a 2 ) 
+ 2C r Hh s (2 (m 2 H (a -5)+ m 2 (a + 6)) - a (a 2 - 5 2 )) , 



G 



K,T(1) 
Hh 



(3.9) 



8C Hh s (2 (s - {m 2 H + m 2 h )) + a 2 - 5 2 ) + WC^a5 (m 2 H - m 2 ) , (3.10) 



G K /h (0) = ± (c^ 2 ((- - *) 2 - 4^) ((a + *) 2 - Am 2 H ) 

- SC^aS (ml - m 2 h ) (s (2 (m 2 H + m 2 h ) - a 2 + S 2 ) 

- 2a6(ml-ml))) , 



G 



K,Y(0) 
Hh 



s((5 2 - 4ml) (S 2 ~ 4m 2 h ) + ^ (™l - m 2 ) 
+ 2a 2 (2 (m 2 H + m 2 h ) + 5 2 ) - 3a 4 ) 
- 4a5 (m 2 H - m 2 h ) (2 (m 2 H + m 2 h ) + a 2 + 5 2 ) , 

and here we have defined the coupling functions as 



(3.11) 



(3.12) 



C 



r 



C 



C 



Hh 

T,l 
Hh 

7,2 
Hh 



U Hh + ° 



Hh ' 



cf xih cf x ° jH (cf x ° ih y 

cf xlH cf x ° ]h (cf x ° iH ^ 



c x ' x ° jH 



(3.13) 
(3.14) 

(3.15) 



C Hh 



c 



Hh 



c- 



y,2 

Hh 



r y,i , r y,2 

U Hh + U Hh ' 



,x?x°A' 



iS 5 IS 1 J \ S /' 



(3.16) 
(3.17) 

(3.18) 
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Higgs (h, H)-neutralino cross term: 



LO 



(h,H- x °) _ 



hH 



=1 r=h,H 



v°v°r 

s — ml + im r T r 



CfG C h + 



ijl n C- 



hH , 



, (3.19) 



= - 4sa+ (4m x o (s 2 - so (a + m^p)) + s (a (2 (m 2 H + m 2 h ) + <5 2 ) - a 2 ) 



G 



C- 
hH 



S(s- a 2 ) (m 2 h - m 2 H ) T s, a, 5, m 2 h ,m 2 H , m 2 



(3.20) 
(3.21) 



where the coupling functions are 



a 



s 



c 



s 



r ijl_ r X?X?H X^h r xU> r ^H 



(3.22) 
(3.23) 



A i Aj 



AA 



Contributions to to come from s-channel Higgs boson exchange, t- and u-channel neutralino 
exchange and cross terms 



~(h,H), r (x°),Jh,H-x?) _ 



AA 



(3.24) 



S-channel CP-even Higgs bosons (h,H): 



~(h,H) _ 1 



LO 



AA 



4 ^ 

r=h,H 



(jrAAfjXiXji 



s — mi + im r T r 



(s - a 2 ) 



(3.25) 



T- and U-channel neutralino: 



*W = \ t ^ A (ct' A ) ' ( Vx? 4 A + ^ + <<) . 

fc,Z=l v ' 

(3.26) 

where 

4l^= (s-a 2 )(r -^o)(s,(T^,mi,mi,m2o,m2o), (3.27) 
= (2^ + | (a 2 - 4m 2 - £ 2 ) T 

+ o-J^i + a (s - a 2 ) ^o) («, cr, 5, m^, m 2 A , m\ , m 2 ) , (3.28) 
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AA 
kl 



(-% - \ (2, + a 2 - Am\ - 5 2 ) T, 

2 _ S 2 ) 2 - 8 m 2 A (a 2 + 5 2 ) + 16ml) (% + 3>o) - 3^2 



16 



(3.29) 



Higgs (h,H)-neutralino cross term: 



_(fc,H-X°) 



i 4 

i5> 



E 

r=h,H 



qAAt (jX-i Xj r 

s — m 2 . + im r T r 



D l> D 



x |^4(7 + (^4m x o - o- 2 + crm x p^ + a 3 - a {km\ + 5 2 )) 
x T s,a,5,m A ,m A ,m 2 



(3.30) 



Y P V ? 

Ai Ai 



hA 



Contributions to d> come from s-channel Z and Higgs boson exchanges, t- and u-channel 
neutralino exchange and cross terms 



■(A) 



W x°x°->^A - u hA + uJ hA + uJ hA + w M 



( z ) ^ ,~.(*°) ^ ,~,(^) ^ ri A -*°) ^ ri 2 ^) . 



/iA 



hA 



(3.31) 



S-channel CP-odd Higgs boson A: 



(A) 
hA 



c hAA c XiXiA 



s - m 2 A + imA^A 



(s + S 2 ) ; 



(3.32) 



S-channel Z boson: 



,~,(Z) _ 1 
" 6 



QhAZ (jXSXjZ 



4 2 



m 7 s 



s — m 2 z + im^r^ 
x |m|s (2s + 5 2 ) (m\ + (s - m 2 h ) 2 - 2m\ (s + m£)) 
~ ( 4( ^ 2 ( m A ~ m 2 h ) 2 m% - m 2 z ({rn 2 A - m 2 h ) 2 m| 
+ 25 2 (3(mi-m 2 ) 2 +(m 2 A + m 2 )m!))s 

+ (6 (mi - m 2 h ) 2 m| - 4 (rr& + m 2 ) m 4 + 5 2 (3 (m\ - m 2 ) 2 + ml)) s 2 
- (3 (mi -m 2 ) 2 - 2ml) s 3 ) a 2 }; 



(3.33) 
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T- and U-channel neutralino: 



J hA ~ V x* Xi *' 4 H feZ 
fc,z=i 



) , (3-34) 



where 



= ^ (« - <5 2 ) (C^4,r r - C+ A Y y ) (s,a,5,mlm%m 2 x o,m 2 x o) , (3.35) 
J h u = (c^ST, + 1 (2C+ >r a (mi - m 2 ) ( s - 5 2 ) 

- (7- A r <5 (a (2 {m\ + m 2 ) + a 2 - 5 2 ) + 2ab (m 2 -mi))) % 

+ ^ C hA,Y a ( m h ~ m A) ( s ~ S 2 ) y<^j (s,v,5,m 2 h ,m 2 A ,m 2 x o,m 2 x o) , (3.36) 

K hA _ ( l r + r r K,T{\) r r K,T(0) T , 1 r + v 
K kl ~ \-2 C hA,T J 2 + G hA 2 ^ +Lj hA J 0+2 L hA,Y y2 

+ J^ C hA,Y aS ( m A ~ m h) yi + G hA {0) y^j ( s > ^ 5 > m l,m 2 A ,m 2 x o,m 2 x ») , 

(3.37) 

and where 

<A T(1) = ( 2C W* 5 K " ml) + C+ AT (-2s 2 + s (2 (mi + m 2 ) + a 2 - S 2 ) 

+ 2a5(m 2 -m 2 A ))) , (3.38) 

G ^ T(0) = 32^ (- 4C, M,T^ ( m A " m^) {s (2 (mi + m 2 ) + a 2 - 5 2 ) 

+ 2ab (m 2 h - mi)) + C+ AT (&aS (m\ (s - aS) -m 4 h {s + aS)) - s 2 (a 2 - 5 2 f 
+ 4m 2 h s (a + 5) (s (a + 5) + cr<5 (5 - a)) 

- 4mi (4m 2 ( s 2 - a 2 5 2 ) - * (s (a - 5) 2 + 5a 3 - <5 3 a) ) ) ) , (3.39) 

G K hA Y{0) = (8a5 ( S (m£ " ™D + ^ K + O) + ^ (- 4 + 2- 2 ^ 2 " 35 4 ) 

+ 4m 2 h s (s (5 2 - 2a5 - a 2 ) + aS (a 2 + 5 2 )) 

+ 4mi (4m 2 (s 2 - a 2 5 2 ) + s (s (5 2 + 2a<5 - a 2 ) - o5 (a 2 + 5 2 )))) , (3.40) 
and the coupling functions are 



hA,T — b S 



Cf (cf ) , (3.41) 
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c, 



± 



hA,Y 




Higgs (A)-Z cross term: 



~{A-Z) u 



(3.42) 



(jhAZ(j x °i x< j Z 



m A + iniA^A ) s — m z + imz^z 



a (s — J 2 ) (s — m|) (m^ — m^) 



Higgs (A)-neutralino cross term: 



(3.43) 



hA 



' (jhAAgxIxlA \ * 
s - m? A + iniA^A j 



hA^A.hA ^hA u A,hA 



(3.44) 



r< c hA _ 

u A,hA — 

r< c hA _ 

u A,hA — 



2s 

1 

4? 



(s — 5 2 ) (^2sm x o + a (m 2 A — m\)^j T s, a, 5, m 2 h ,m 2 A , m 2 (3.45) 



5 1 4s + I s (4m o - 2 (m 2 A 



+ 5 2 -a< 



+ 2a5 (m A - m 2 h )) T s, a, 5,m 2 h ,m A ,m 2 x o ) 
Z-neutralino cross term: 



(3.46) 



(z-x°) 



hA 



r< c hA 

^Z,hA 



l=i 



Ch AZ C X°X°Z y , 

- m 2 z + imzTz ) \ C ^z, h A + ^hA^z,hA 



(3.47) 



+ 



Y s (s (2 (s-m 2 A - m 2 h ) + 4mJ„ + c5 2 ) + 2ab (m\ - m 2 h ) - a 2 s^j 

-«5 4 ) 



8m\? ( 16 K? " ( m x? " m h) + 16m lf 

— 8as (m 2 A — m 2 ) ( m x? ( m z ~~ s ) ( s ~ ^) 
+ m 2 z 5 (m 2 A + m\- s) - 2m 2 om|5) 

+ 4a 2 (s 2 ((mi - m 2 ) 2 - m| (mi + m 2 + 2m 2 ( ,)) 

— (m 2 4 -m 2 ) 2 (,-2m|) ( 5 2 ) 

— 4m 2 z sa 3 5 {m 2 A — m^) + m 2 z s 2 a 4 ) T s,a,5,m\,m 2 A: m\ 



(3.48) 
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{~i C hA 

u Z,hA 



\-5 [a [m\ - m 2 h ) (s - m|) - 2m x om|s) 
\-^5 (2m x om 2 z s 2 - 2 {m\ + m 2 h - s) - 5 2 ) 



m z s 



4m z s 



+ as (771/4 ~~ m h) [2 y~^ m2 z ( m A + m h ~ ^ m x° 
+ s (jn A + 777^ + m 2 z — 2777^,0^ + Am x om 2 z 5 — s5 2 ^j 
— 2a 2 (^s 2 m x om 2 z — (m\ — ml) 2 (2m 2 z — s) 5*J 
+ sa 3 (m 2 A - m 2 h ) (s - 2m|)) T s, a, 5, m 2 h ,m 2 A , m 2 ^ 

where the coupling functions are: 



(3.49) 



± r x°^h x o x o A x o jX o A o x o h 



(3.50) 



Contributions to the HA final state use the same expressions as above, but with mh 
replaced with 777/7 and h replaced with H in the couplings. 



Contributions to u) come from s-channel Z and Higgs boson exchanges, t- and u-channel 
chargino exchange and cross terms 



_ .(MO , _(Z) , ~(x±) (MZ-x*) , ~(z-*±) 

^X^X^H+H- - ^H+H- + W H+/7- + + W //+/7- + 

S-channel CP-even Higgs boson (h,H): 



(3.51) 



~{h,H) 
H+H- 



r=h,H 



C 



rH+H- 



CXiXjr 



s — m 2 . + im r T r 



(s-a 2 ); 



S-channel Z boson: 



(Z) 

H+H- 



C 



ZH+H- (jXiXjZ 



A 



s — 777,7 + imzfz 



(s-Am 2 H± )(2s + 5 2 ) (s-a 2 ) 



T- and U-channel chargino: 



(3.52) 



(3.53) 



k,l=l 



Xk 



(3.54) 



where 
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+ (s-6*)-Cf + l -_ (s-a 2 )) 

x Vo) (s,cr,5,m 2 H± ,m 2 H± ,m 2 x ±,m 2 x ±) , (3.55) 

jH+H- _ l (( ( Mij}kl-+ Mij}kl+-\ 
J kl ~ 2\\\ H+H ~ +L, H+H- ) a 

- (S^t" " Df: i H t-) 5 (4m 2 H± + a 2 - 5 2 )) T 

+ (^--^;-)*(-^)) 

x y ) (s,a,5,m 2 H± ,m 2 H± ,m 2 ± ,m 2 ± ) , (3.56) 

- 2 (z)<^_ + + ( S - 2m!*)) * 

" 16 V H + H- +L/ H+H- ) 

x (l6m^± - 8m 2 H± {a 2 + <5 2 ) + (a 2 - ,5 2 ) 2 ) T 

+ ^ ((^S- + ( 8 -h ± (* 2 - - 2 ) + 2 (- 4 - * 4 )) 

+ (^S- " (I6m| ± + 25V - a 4 - <5 4 )) 

x 3>o) (s,cr,S,m 2 H± ,m 2 H± ,m 2 ± ,m 2 ± ); 

(3.57) 



Higgs (h,H)-chargino cross term: 



- (h,H- 



i 2 



rH+H- 



fe=l 



s — to 2 . + im r r r 
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X 

x T 



{ c T+h- (- 4cj + a ( 4m H± - 4 ™'± + ^ 



s,a, 5, m H ±,m H ±,m ± 



C l ^ H ^4m x ± (s-a 2 )^ s,a,S, m 2 H ± , m 2 H± , m^± 



(3.58) 



Z-chargino cross term: 



H+H- 



5> 



k=i 



(jZH+H- (jXiXjZ \ 



s — m z + imz^z 
^(2s-4m 2 H±+ 4 m 2 xi -a 2 -S 2 ) 



+ ^(a-o 2 ) {s-4m%±) 



+ 



4m 2 H± +4m 2 ^ -oZ-py^j 



s,a,5, m H± ,m H± ,m ± 



where the coupling functions are 



(3.59) 



^djkl± 
U H+H- 
f=,{ij}kl± 

fiijkl-\ — 

f-iijkl — h 
^H+H- 
^{ij}kl-\ — 
^H+H- 
-~,{ij}kl- + 



c 



^H+H- 
^H+H- 



f-djl,± 
^H+H- 

s~<jil,± 



C 



ijk,+ f-djl- 
H+H-^H+H 



(~iijk,+ 
^H+H- 



H+H- 



C 



ijk- 
H+H 



' H+H~ 



, s~ijik,± 
+ ^H+H- 

, r ijik,± 
+ U H+H- 

, n jik,+ 
+ ^H+H- 
r jik,- 

, n jik,+ 

n iik- 



^H+H- 
^H+H- 



■jil- 
H+H- 

-^H+H- 

-^H+H- 



C 



H+H- 



(3.60) 
(3.61) 
(3.62) 
(3.63) 
(3.64) 
(3.65) 



The same expressions are used to define the D couplings, but with the replacement 
C —>■ D. We also need 



D 
D 



[ij]kl+- _ n ijk,+ n jil- 
H+H- ^H+H-^H+H- 
\ji]kl+- _ n jik,+ pfijl- 
H+H- ^H+H-^H+H- 



^H+H-^H+H- 
njik,- n ijl,+ 
' ^H+H-^H+H- 



(3.66) 
(3.67) 



All of these souplings are based on the combinations 



C 



ijk,± 
H+H- 



C X°xtH-\ c X$xtH- ± c xMH-\ c X { }x + k H 



n ijk,± _ ( r xUi H ~ 

U H+H- ~ °P 



°5 



X jXtH- 



C xUiH-\* c X°xtH- 



(3.68) 
(3.69) 
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oo vW + H 



Contributions to uj come from s-channel Higgs boson exchanges (both CP-odd and -even) , 
t- and u-channel chargino exchange and cross terms 

~ (A) . ~(h,H) . -(**) . -{A-X^ . -{KH-x*) fonn\ 

u x o x o W+H - = ^ W > +H _ + w^'+jy- + ^ W+ ' H _ + u y w+H J + uj v w+h _ : (3.70) 



S-channel CP-odd Higgs boson A: 



~ (A) 



C 



AH-W+ (jXiXjA 



s - m A + ituaTa 



(s - 5 2 ) (s 2 - 2 (m 2 H± +m 2 w )s+ ( 



m 2 H± - m 2 w 



2m 2 w 



(3.71) 



S-channel CP-even Higgs boson (h,H): 



U W+H~ ~ 1^ 
r=h,H 



rirH-w+ rtXiXji 



s — m 2 . + im r F r 



2m^ 



2 2 \ 2 



(3.72) 



T- and U-channel chargino: 



where 



! V \m rr, T WH _J_ ™ _l_ l/TO 

— 2^ |"\: "V '/./ • "'. •'/•' + K ' 



~Xk 



'kl 



where 



t wh _ 
1 ki — 



r I,T(2) r r I,T{l) r r I,T(0) r 



(3.73) 



J kl — 



- G I ^ 2) y 2 + G I ^ 1) y l + G I ^ ) y Q ) (s,a,5,m 2 H± ,m 2 w ,m 2 x ±,m 2 xt ) ,(3.74) 

( r J,T{2) r r J,T(l) T , r J,T(0) r 

+ G J // ) y 2 + + GiKJ 0) ^o) (5, a, <5, < ± , m 2 w , m 2 ± , m 2 ± ) 

^Wif _ ( r K,T{2) r r K,T{l) r r K,T(0) T r K,Y(2)„ 

+ Gwl^yi + GwH {0) y<>) (s,*,5,m 2 H± ,m 2 w ,m x ±,m x ±), (3.76) 



1,(3.75) 



G 



I,T(2) 
WH 



D 



c,+ 

WH,ij 



D C,+ \ 



-2a + (6 - a) 2 ) 



(3.77) 
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+ 2 (m 2 w (s - 2o5) + m 2 H± (s + 2<r<5)))) , (3.78) 

G^f = ^ (s 2 D C w + H, i3 (~^w (8* + ^ - 2a5 - 95 2 ) 

+ Am 2 H± [Am 2 w + (a + <5) 2 ) + (a + 5) 2 (-4s + (a + <5) 2 ) ) 

" rfj&ji (« 2 ("4- + (S ~ <??) (5 - a) 2 
+ 16m w aS (—s + a 5) + 16m H ±aS (s + a 5) 

- Am 2 w s (8s 2 + 2a5 (5 - a) 2 + s (3a 2 - 2ab - 95 2 )) 

+ Am 2 H± [Am 2 w (s 2 - 2a 2 5 2 ) + s (s (a 2 - 6ad + 5 2 ) + 2ad (a - 5) 2 ))) 

- A8m 2 w s 2 (a 2 - 5 2 ) ((D+ H . jk )* D WHijl + (D WH , jik ) * D+ HdU ) ) , (3.79) 

G vS 2) = - 2 ~ 2jD W H ff,{ii}) ' ( 3 - 80 ) 

GfeT = " f s { C WH M ~ 2D wh M ) {" + ™w ~ m 2 H± ) , (3.81) 

G ^ 0) = i (( 2D ?&,tt> " C wtm) ("8m^ (* " + 8^ ± aJ (, + a<5) 

+ 4m#± (s + aS) (s [Am 2 w + (a - J) 2 ) - Am 2 w ob} 

+ s 2 ( <T 4 + 6 < 7 2 <5 2 + <5 4 -4 S (<7 2 + (y 2 )) 

- 4771^5 (8s 2 - s (3a 2 + 2a5 + 35 2 ) + aS (a - <5) 2 )) 

+ 2AC^ m m 2 w s 2 (a 2 - 5 2 )) , (3.82) 



\ (° Re ( C WH, {ij} ~ D WH, {ij} ) ~ 6Re ( C WH M ~ D wh M )) ' ( 3 - 83 ) 
^ { Re { C W + H M - D w + h M ) + (s (S-a) + 2m 2 H± a) 

Re { C w + H,m - D w + H,m) + ( s ~ s ) + 2m H^) 

8m 2 w sRe (s C W + HM + a D w + H {ij} ) - Am 2 w sRe (a C w + H {ij} + 5 D w + H [ij] ) 
2m 2 w Re [a (c"+ M - D W + HM ) 
+ * { D WH,m ~ C w + H, m ) )<r6{* + S)), (3.84) 

G ^ 0) = 32^ ( Re ( C WH,m ~ D w + H,m) ° ( 8m w (-4s 2 + s5 (a + 5) + ob 2 (a + 8)) 

- Am 2 w (a + 5) (s (S - a) (s + S (S - a)) + 2m 2 H± 5 (s + 2a8)) 

+ (a + 5) (Am 2 H± s5 (a - 5) 2 + 8m A H± 8 2 a + s 2 (a - 5) (a 2 + 35 2 ))) 

+ 12i?e (C~ + H + D w + H m ) m 2 w sa (2m 2 H± (2s + 6 (a - 5)) 

- (a - S) (s (a + 5) + 2m 2 w 6)) 

~ Re { C w + H M - d w + h M ) 6 ( a + 5 ) { sa ( 2m ^± + (°- " ^) 



r 1 



J,T{2) 
WH 

J,T(1) 
WH 
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+ Am 2 H± a 2 8 + s 2 (a - 8)) 

+ 8m^ (-4s 2 + sa(a + 8) + a 2 8 (a + <5)) 

+ (a + 8) (Am 2 H± sa (a - 5) 2 + 8m 4 H± a 2 8 + s 2 (5 - a) (3a 2 + 5 2 ))) 

+ 12Re (C~ + HM + D w + H [ij] ^j m 2 w s5 (2m 2 H± (2s + a (8 - a)) 

+ (a -5) (s {a + 5)+ 2m 2 w a))) , (3.85) 

<®P = \ (* + D^'l } ) + S ((Cjf+H + D^)) ' (3 - 86) 

+ a8 (s + m 2 w - m 2 H± ) [a + 8 D^~}^ 

G ^ 0) = 64^ (a (5fi2 + ^ (CT " 55) " 52(j2) 

+ 8m 4 H± (s 2 + s5{a-5)+ 5 2 a 2 ) 

+ Arn^s (8 - a) (8a (a - 8) + s (3cr + (5)) 

+ s 2 ( 8s 2 _ As ( a 2 + 3^) + ^4 + ^2^2 + ^ 

+ 4m^± (4m^ (3s 2 - 2s<5 2 - a 2 5 2 ) + s (-4s 2 

+ (5 (s 2 (8 (m 4 H ± + m\y - 2s (m^± + 3^^) + s 2 ) 

+ 45 2 (m^± + 9m 2 w - s) + 8 4 ) 

+ 4sa8 (m 2 H± - m 2 ^) (2 (m 2 H± + m 2 ^ - s) + <5 2 ) 

+ 2d 2 (2s (-3s 2 + s (5m|, ± + m 2 w ) - 2 (m 2 H± - m 2 w f^ 

+ 5 2 (3s 2 - 4s (m 2 H± - m 2 w ) + 4 (m^± - m 2 ^) 2 )) + sV 

+ a (s (8s (m#± + m%[ - 2s (m 2 H ± + Sm^) + s 2 ) 

- 4<5 2 (2 (m 2 H± - m 2 w f - s (5m 2 H± + m 2 w ) + 3s 2 ) + s<5 4 ) 

+ 4scr<5 (m 2 H ± - mly) (2 (m 2 H ± + m 2 ^ - s) + 8 2 ) 

+ 2<r 2 (2s 2 (m^± + Vmly - s) 

+ 5 2 (3s 2 - 4s (m^± - m 2 w ) + 4 (m^± - m 2 ^) 2 )) 

+ 4s<x 3 5 (m 2 ^ - m 2 w ) + s 2 a 4 ) D#+-> } 

+ 8 (8m 4 H± (s 2 + sa(8-a) + a 2 8 2 ) - 8m^ (5s 2 + sa (5 - 5a) - a 2 S 2 ) 

- Amlys (8 - a) (s (a + 38) + a8 (8 - a)) 

+ s 2 (8s 2 - 4s (3a 2 + 8 2 ) + a 4 + 6cr 2 J 2 + 8 4 ) 

+ 4m 2 Y ± (4771^ (3s 2 - a 2 (2s + J 2 )) 
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+ s (_ 4s 2 + s (5a 2 - 2a5 + 5 2 ) + a5 (5 - a) 2 ))) D^+J^ , (3.88) 
Gwh 2) = ~l {{ D WH,i, + Vfylu) + + ^ + ^ 2 ) 

G ^ (1) = lb (^« a ( 16m ^ + (" 2S + (5 " a)2 ) (<5 + CT)2 

- 2m 2 w (8s + (55 - 7a) (5 + a)) + 2m 2 H± (sm 2 w -(5 + a) 2 )) 

+ D°f Hjji (l6m 4 w (s - 25a) + s (5 - a) 2 (-2s + (5 + a) 2 ) 

- 2m^ (8s 2 + 25a (5 + a) 2 + s (55 2 - 65a - 7a 2 )) 

+ 2m 2 H± (25a (5 + a) 2 + 8™^ ( s + 25<r) - s (5 2 + 6<fo + a 2 ))) 

+ 2Am 2 w s (5 2 - a 2 ) ( (D WHijk ) * D WH ijl + D WH,jii)) > (3-90) 

- 16m^ (2s 2 (5 - a) 2 - s<5a (5J - 7a) {5 + a) - 5 2 a 2 (5 + af^j 
+ 16m A H± 5a ((5 + cr) 2 (-s + <5cj) + 8771^ (s + <5cr)) 

+ Am 2 w s {5 - a) 2 (-25a {5 + a) 2 + s (5 2 + 65a + <r 2 )) 
+ 4m 2 H ± (-s (s - 25a) (5 2 - a 2 f + 32m^ (s 2 - 25V) 

- Am 2 w (asSo {5 - 2a) (5 + a) + 25 2 a 2 (5 + a) 2 + s 2 (<5 2 + 6<fo + ° 2 )))) 
+ D WH,ij s2 ( Am H± - ( S - °?) ("32m^ + Am 2 w {5 + a) 2 + (5 + af) 

- 192 (^Hjifc)* ^W,^^ Kr± - - ° 2 )) . (3-91) 

« (2) = 5 (-*4<h >W } + C£S, W , (4- - - 2 + * 2 ) + ZD?- >{y} (a 2 - 5 2 )) , 

(3.92) 

= — {- 8m W C WH,[ij} + ( 2 - C ' 2 yif,[ij] ~ C WH,[ij]) S ~ 6<2 ) 

+ m 2 w (-4sC%- M + (a 2 - 5 2 ) (2D^ HM - 

+ < ± (Sm^^- ^ - (a 2 - 5 2 ) (2D^ HM - C#£ >[y] ))) , (3.93) 

G WH (0) = Y2^2 ( fi2 { l2Sm H± m W C WH,{ij} ( S - 

- 16 (-m|, ± C^- + m 2 H± (sC°'- >{ij} + (7C£'- {ij} - 4D& Hm ) 
+ m 2 ^ (m 2 ^ (c^J H{ij} - 6D^ H{ij} ^ + s (c^ 

+ 2 ^, { y } ) ) ) ) ^ + 454 K (<"..;;, " GD WH, {ij} ) 
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+ m 2 H± (c$ H){ij y 2L >WH,{ij}) + 2sD WH,{ij}) 
+ 56 ( C WH,{ij}- 2D WH,{ij})) 

+ 4scj5 (m 2 H± - m 2 w ) (2 (s - m 2 ^ - m^±) - <5 2 ) (8m^C^'~ { . j} 

+ ^k lW > " C WH,m ( 4 * + ^)) + ^ (!6« 2 (C%~ m « ± + 5m 4 w 
- 3m^s - m^± (s + 3m^)) + 2m 2 v D 1 f H{ij ^ (s - 3m^ - 2m#±)) 



+ 4m^s (m^ - m#±) - 4s 2 (m 2 ^ + m#±) + s 3 )) 

+ 54 " 2Z> w^,{ii}) ( 5fi2 + 8s (™w ~ m 2 H ±) + 8 (m 2 w - m 2 H± f) 

8sa 3 5 {m 2 w - m 2 H± ) (-C^ { . j} (3s - ^m 2 w - m 2 H± ) 



+ 



+ 2D WH,{ij} ( S ~ m W ~ m H± 

+ a 4 (As 2 



-m 2 H± (c^ H {ij} 2D 2 ^ H {ij ^ + 2s (c^ H {ij} DwH,{ij}) 

~ 52 { C WH,{ij} ~ 2D WH,{ij}) ( 5g2 + 8s ( m w ~ m H±) + 8 ( m w ~ mjj±)' 
+ (_ S V + AsoH (m 2 w - m 2 H± )) (c%~ m - 2 D 2 W ' H ^ } )) , (3.94) 
where we have used the following combinations of gauge functions: 

C WH,{ij} = C WH,{ij} + C WH,{ij} ' ( 3 - 95 ) 

C WH,[ij] = C WH,{ij} ~ C WH,{ij} ' (3.96) 

D WH,{ij} = D WH,{ij} + D WH,{ij} ' (3.97) 

D WH,[ij] = D WH,{ij} ~ D WH,{ij}> ( 3 - 98 ) 

^WH,{ij} = {pWHjik) ^WH,ijl + {pWH,ijk) 

^WH,[ij] = {pWH,jik) ^WH,ijl ~ {pWH,ijk 

C]VH,{ij} = (f'WHJik) ^WH,ijl + [pWH,ijk 

^WH,[ij] = ipWHjik) ^WH,ijl ~ {pWH,ijk) ^WHJil ' 

D WH,{ij} = (y D WH,jik) D WH,ijl + iPwH,ijk) 

D WH,[ij] = iPwH,jik) D WH,ijl ~ ( D WH,ijk) 



D WH,{ij} ~ ( D WH,jik) D WH,ijl + ( D WH,ijk) 
D WH,[ij] = ( D WH,jik) D WH,ijl ~ ( D 



'WH,ijk 



CwHjil > 


(3.99) 




(3.100) 




(3.101) 




(3.102) 


D WH,jil ' 


(3.103) 


D WH,jil ) 


(3.104) 


D + 


(3.105) 


D + 


(3.106) 
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^WH,{ij} 
^WH,[ij] 
^WH,{ij} 



a 



WH,[ij] 



D 



WH,{ij} 



D 



WH,{ij} - 



^WHJik J ^WHJil + \pWH,ijk ) ^WH,ijl 

{pWH,jik) ^WHJil ~ {pWH,ijk) ^WH,ijl 

(pWH,jik) ^WH,jil + (pWH,ijk) ^WH,ijl 

CwH,jik) ^WH,jil ~ (pWH,ijk) C WH ijl 

\\ll,jik) D WH,jil+ {j^WH,ijk) D WH,ijl 



( D WH,jik) 



D 



WHJil 



( D WH,ijk) 



D 



WH,ijl 



DwHJik) DwHjil + \^WH,ijk) ^WH,ijl 



( D WH,jik) 



D 



WH,jil 



( D WH,ijk) 



D 



WH.ijV 



D 



2,± 

WH,{ij} 



D ± ) D ± + ( D ± ) D ± 

^WHjik) ^WH,ijl^ X^WH^jk) ^WHJil 



D 2,± _ ( D ± y D ± _( D ± y D 



± 

WH,jil 



a 



D,± 



WH,{ij} 



a 



D,± 
WH,[ij] 



p± , p± 

1 ijkl ~ 1 jikl ' 

p± _p± 



ijkl 



D WH,ji — ipWHjik) ^WH,jil + {pWHjik) D 

{ D WH,jik) 



P ± ~ 



(CwH,jik) ^WH,ijl + 



± 

WHjil 
± 

WH,ijV 



U WHjik - U 5 1 A 



±c 



xUt^ ( c xUl w± 



±c 



v 



xtx°H± ( c xUl W± 



V 



Higgs (A)-chargino cross term: 



W+H- 



i=i 



(jAH-\¥+ (jXiXjA ' 

s — m? A + iniA^A 



1m 2 



Di 



n H- ^WH,^} 
U WH,{ijY J WH 



I n i<.- r D WH,M , n H,+ r D WH,{ij] n H,+ r D WH,[ij\ 
+ U WH\ijY J WH + U WH,{ijY*WH U WH\ijY J WH 



, (3.123) 



where 
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Wl 



D 

WH,{ij} 



WH,[»j] 

WH 



-m x +cj (s - <5 2 ) (( s - TO ^) 2 + m H± 
2m 2 H ± (s + mw))? 7 s,a,5,m 2 H± ,m 2 v ,m 2 + 

L A; 

^m x +<5(s + m^-m^±) 
4s + (4m 2 + s + s (2 (s - m 2 ^ - m^±) - 5 2 ) 



WH 



— (4s (2s {s-mly- m 2 H ±) - S 2 (s + 



4s 



m^±)) 



- (s (-8s (2m 2 v m 2 H ± + m 2 + (s - m 2 ^ - m 2 H± )^j 
+ 25 2 ((m 2 ^ - m 2 H± ) (2m 2 + + m^± - 3m 2 w ^ 



<5 4 (s + m 2 ^ - m 2 H ± 



+ s ( 2m 2 + + m%± + Sm 2 ^ 

+ 2cT<5 (m 2 ^ - m 2 H± ) (2s (m#± + m 2 ^ - s) 

+ 5 2 (s + m 2 ^ - + mly)) 

- a 2 (2s {{m 2 H ± - mly) 2 - s (m 2 H ± + m 2 ^) 

- 6 2 (2 (m 2 H± - m 2 w f - {3m 2 H± + 5m 2 w ) + s 2 ))) 



D ~ 



x T 

1 



s,a, 5, m 2 H ±,mw, 



m 



(3.124) 



2a5 (m 2 H ± — m^y) — so -2 ) T s,a,5,m 2 H ±,m 2 v ,rn 2 + J , (3.125) 

L " ' X; J / 



(3.126) 



—(j(5 ^-4m^s + Amlys 2 + 4m 2 + s (s + m 2 ^ - m^±) + 2m 4 y (5 2 

— hrn^b 2 + s 2 5 2 + 2m 2 ycr<5 (s + m 2 ^) — sct 2 (s + m 2 ^) + 2m 4 ?± <5 (a + 5) 
+ (s (a - 35) (a + 5) + 4™^ (s - <5 (a + 5)))) ; (3.127) 

Higgs (h,H)-chargino cross term: 



where 



_(ft,H- X ±) 



Z=l 



Is — m 2 . + im r T r 



m 



w 



^WH,{ij} U WH ^WH\ijY^W 



WH,[ij] 



~ r ^WHJii^WH ^WH 



'WH\ijY JX WH 



WH,[ij] WH 



(3.128) 



— m x ±(J ( s + m VK ~~ m H ± ) + \^ m x ±,S ~~ ^m^/S + 2s — S(5 
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r H,+ 



,11, 
'WH 



+ 2my/(j5 — sa 2 — 2m 2 H± (s + <7<5)) T s,a,5,m 2 H ±,m 2 v ,rn 2 ± J , (3.129) 

L X; J / 



— m x ±5 (2s - a (a + 5)) (mj I± + (s - rn^f 

- 2m 2 H± (s + m 2 v ))j :r s,a,5,m 2 H ±,m 2 v ,rn 2 ± 

L M . 

- ( 4s ( 8fi2 ( m ///± + m H / ~ S ) ~ 2s0 ~^ ( 3m ///± + m W ~ 3s ) 

+ (T 2 {s + m 2 w - m 2 H± ) (4s - 3<5 2 )) 

+ ^8s 2 ^— 8rn 2 H± m 2 v s + 4m^±s (wi^i + ra 2 ^ — s) 

<5 2 (s (m#± + mfy) - (m#± - m 2 ^) 2 )) 

- 2sct(5 ( 4m 2 ±s (3m 2 H ± + m^y — 3s) 

- 4s (-2m^± + mw - ntfys + m 2 H± [hm^r + 2s)) 
+ 5 2 ( s 2 + s (Sm 2 ^ + m 2 H ±) - 2 (m^± - m^f 



+ S(7 

+ s 



+ 



^8s ((m 2 ^ - m 2 H ±) ^2m 2 ± + - 3m^^ 
2m 2 ± + m 2 H± + 3777^ ) 



+ 45 \ {m H± -m w ) \^ m x ± +5m H ± 



4m 2 w 



- s (srn 2 ^ ± + 6m 2 H± + Sm 2 ^^ + s 2 ^ + 3<5 4 (s + rriiy - m 2 H± )^j 
+ 2o- 3 5 (s (8 (m 2 H± - m 2 w f - 3s (3t77^± + m 2 w ) + s 2 ) 

- 35 2 ^2 (m 2 H± - m^f - s (3m 2 H ± + mfy) + s 2 ^ 

+ sa 4 (m 2 H ± - - s) (4s - 35 2 )) J 7 s, a, 5, m 2 H ± , m^, m 2 ± ) , 

L M J / 



T a6 (4s (2s - (7(5) (s + 



+ (2s 2 (4777^ (s - m 2 ^ + 777^ ± ) - 5 2 (s + m 2 ^ - m^±)) 

+ S(7<5 (4777^± + 8777^ + 5 2 (s + 777^) - 777^± (4s + 12777^ + 5 2 )) 

+ 2d 2 (s (s 2 - s (3777^ + 5777^) + 2 (m 2 ^ - m 2 ^) 2 ) 

- 5 2 ^s 2 - s (3m 2 H ± + m 2 ^) + 2 (m 2 H ± - m 2 ^) 2 ^ + s<t 3 5 (s + m 2 ^ 

+ 4777 2 ±s (2s - (7(5) (s + - m 2 H ±)^j T s,ct, 5, m#±,m^,m£± ) 

where the coupling functions are: 



(3.130) 



(3.131) 



m 2 H± ) 
(3.132) 



D 



H- 

WH,[ij] 



D 



WH,{ij} 



D WH,ijl + D WH,jil 
™WH,ijl ~ ^WH,jil 
D WH,ijl + D WH,jil 



(3.133) 
(3.134) 
(3.135) 
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D WH,[ij] = D WH,ijl ~ D WH,jil ' (3.136) 

^WH,{ij} = C WHijl + C WHjU , (3.137) 

< ~'WH,[ij} = C W H,iji~ CwHjiV (3.138) 

^WH,{ij} = ^WH,ijl + ^WH,jil ' (3.139) 

^WH,[ij] = ^WH,ijl ~ CwH,jil-> (3.140) 



= Cf^"- (c^-)'±C^>«- (cf^J , (3.141) 

= C***- (cf^- ) • ± of*"' (c^ W - )". (3. 142) 

Contributions to uj come from s-channel Z and Higgs boson exchanges, t- and u-channel 
neutralino exchange and cross terms 

~{A) . ~(Z) . . ~{A-Z) . ~{A- X °) . ~{Z-X°) (oaao\ 

^xh^zh = "zh + ^zh +"zh + w k + w k + w k • ( 3 - 143 ) 
S-channel CP-odd Higgs boson A: 



(A) 



s - mi + iniA^A 



(s - 5 2 ) (s 2 - 2 (m| + m£) s + (m? 



m 



a- 



2m| 



(3.144) 



S-channel Z boson: 



~(2) 
u Zh 



QiZZhfjXiXj- 



s — m z + imz^z 



1 



2s m z [m\ + m% + 10s m| + s 2 — 2m^ (s + m 2 ^)) 



24m|s 

+ cr 2 (3s 4 - 12s 3 m| + 16s 2 m| - 32s m% + m| 
+ (m^ - 2m 2 t (s + m|)) (m| - 6sm| + 3s 2 )) 

+ — (l2sVm| - 3s V 2 + s 2 m| (s - 19cr 2 ) + m% (s - 4a 2 ) + 2s m% (5s + a 2 ) 
+ (m£ - 2m 2 h (s + m|)) (m| (s - 4a 2 ) + 6s a 2 m 2 z - 3sV))) ; (3.145) 
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T- and U-channel neutralino: 



4 

4? = A- E hx^x?4? + + *« fc ] > ( 3 - 146 ) 



2m z fc 7=i 



where 



4? = (- (4$' + d«>») r 2 + g^t; + 4£ (( H + (4& H + 4& H ) ^ 

- Y s {c^ + C^)aS(s + ml-ml) yi 



+ 4r (0) ^o) (s,cr,S,ml,ml,m^o,rn^o) , (3.147) 
j h - [-[[^zhA + u zhy )° + { 1J zh,A + D zhy) d ) I ^ + L 'zh ^ 

r J,T(0) T ( r {ij}kl r {ij}kl\ v r J,Y(l) v 
+ G Zh J + [yzh,A + C Zh,V ) + t* zh ^1 

+ ^r (0) %) (s,a,S,mlm 2 z ,m 2 x o,m 2 x o), (3.148) 

^ _ ( r K,T{2) r r K,T(l) T , r K,T(0) r 

+ G^f (2) ^2 + G^f + G£f ( °Vo) (s, a, (5, m|, m 2 o , m 2 o), (3.149) 



^ W = £ ( {D»j$ + , (2 (m h 2 + m z 2 - s) + S 2 ) 

+ 2 ((^» + D™* + 4«>» + 4^ } v ) " + 

" (41a + 4^) «) ^ + (41?" + D Wy) *<») . (3-150) 

+ D§^ 5 (m h 2 - m z 2 - s) a 2 (2m z 2 (s-Sa) 

+ 2m h 2 {s + 5a) + s (-2 s + 5 2 + a 2 ))) 

+ D Zh,7 (8m z 4 5a (s + 5 a) + 8m ft 4 5 a (s + 5 a) 

+ 4m A 2 (s + 5a) (s (<5 - ct) 2 + 4m z 2 (s-<5cr)) 

+ s 2 (5 4 + 6<5 2 a 2 + a 4 -4 S (<5 2 + a 2 )) 

- 4m z 2 s (8s 2 + <5(5-fj) 2 a-s (3<5 2 + 25a + 3<i 2 )))) , (3.151) 

4?^= I ^(24C ZM m z 2 S 2 H 2 + . 2 ) 

+ C Z h,3 (8m z 4 5a {-s + 5 a) + 8m h 4 5 a {s + Sa) 

+ Am h 2 (s + 5a) [s (5 - a) 2 + 4m z 2 (s-<5cr)) 

+ s 2 (5* + Q5 2 a 2 + a 4 -4s {5 2 + a 2 )) 

- 4m z 2 s (8s 2 + 5{5-a) 2 a-s (35 2 + 2Sa + 3a 2 )^ , (3.152) 
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G J z T h W = j- g (p% h>1 (2 m h 2 5 a + m z 2 (4 s - 2 5 a) + s (5 2 + a 2 )) 

+ D s z \ 3 (2m h 2 5a-2m z 2 (4s + <5a) + s (<5 2 + <r 2 ))) , (3.153) 

G J z T k (0) = ~ ^ ("24 m z 2 (d 6 z \ 2 (-m h 2 + m z 2 ) + D 5 Z \ A s) (5 - a) (5 + a) 
+ £>zfc,i (~8m z 4 (4s - 3<5 2 + 25<r - 3d 2 ) + s (5 2 - a 2 ) 2 
+ 4m h 2 (5o{5 + <j) 2 + 2m z 2 (6s - 35 2 + 25a - 3ct 2 )) 

- 4m z 2 (2s (<5 2 + a 2 ) + 5a (5 2 - 65a + a 2 ))) 
+ rfzhfi ("32 m z 4 (s - 5 a) + s (5 2 - a 2 ) 2 

+ Am h 2 (5 a (5 + a) 2 -Am z 2 (3s + 25<r)) 

+ m z 2 (40 s (5 2 + a 2 ) -45 a (5 2 + 18 5 a + a 2 )))) , (3.154) 
G? h W = \{{c z f} l {- mh 2 -, mz 2 + s) 

+ C z fv i~ m h 2 + m z 2 + s))a) , (3.155) 
G J z l i0) = ^ (a (-24 C ZM m z 2 (m h 2 (s - 5 2 ) + m z 2 (s + 5 2 ) 

+ (s-2<5 2 ) (s-a 2 )) +C Zht3 {8m h 4 (s + 5a) 

- 8m z 4 (2s + 5 (-35 + a)) + s (8 s 2 + (5 2 - a 2 ) 2 - 4 s (<5 2 + <r 2 )) 

- 4m z 2 (6s 2 + <5a (<5 2 + 6 5 a + a 2 ) -s (95 2 + 25a + 3ct 2 )) 

+ 4m ft 2 (6m z 2 (s-<5 2 ) - (s + <5ct) (4s- (5 + a) 2 ))))) , (3.156) 

<£f (2) = ^ (" ((D%yf + s (8^ 2 - 4s + * 2 )) 

+ 6 (4fi + 4^)s^ 

- (4t? + 4tv) («-45 2 )^), (3.157) 
G«™= ^-(2 s (-6D Zh ^ mz 2 s (5 2 -a 2 ) 

- D 5 z \ 5 a 2 (8m z 4 + m z 2 (6 s - 5 2 - 6 5 a - a 2 ) 

+ 2s (s + 5 2 + a 2 ) +m h 2 (-8m z 2 + 2s + 5 2 + 65a + cr 2 ))) 
+ £>S h>7 (I6m z 4 s (s-<5a) 

+ 2m A 2 5 {5 - a) 2 a - 45 3 a 3 + 8m z 2 s (s + 5 a) - s 2 (5 2 + 45a + a 2 )^ 

+ 2ra z 2 (-8s 3 + 45 3 a 3 -s(5(j (5 2 -W5a + a 2 ) + s 2 (c5 2 + 4 ,5 a + <r 2 )) 
+ s (-2s 2 ( < 5 2 + C T 2 ) -4<5 2 ct 2 (5 2 + cr 2 ) +s (<5 4 + 6,5 2 a 2 + a 4 )))) , (3.158) 

G z ? (0) = ^ fe, 7 (64 mz 6 <5 a (s - 6 a) 
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+ 8m h 4 5a (s + 5 a) (-8m z 2 + (5 + a) 2 ) 

_ s (5 2 -a 2 ) 2 {-A5 2 a 2 + s {5 2 + a 2 )) 

+ 8m z 4 (4s 2 (5-a) 2 + s5a (S 2 - U5 a + a 2 ) + 5 2 a 2 (S 2 + 10 5 a + a 2 )) 

- 4m z 2 (s 2 (5-a) 2 {5 2 + 45 a + a 2 ) + 45 3 a 3 {5 2 + 65a + a 2 ) 

- s5a {5 4 + 125 3 a + 65 2 a 2 + 125 a 3 + a 4 )) - 4m h 2 {32m z 4 (s 2 - 5 2 a 2 ) 

- {S + af (4S 3 a 3 -s5a(5 + a) 2 + s 2 (S 2 + a 2 )) 

- 4m z 2 {s 2 {5 2 + A5a + a 2 ) - s5a {5 2 + 65a + a 2 ) 

- 5 2 a 2 (5 2 + 65a + a 2 )))) + 2a (24Z^ M m z 2 (m ft 2 - m z 2 ) s5 (5 2 - a 2 ) 
+ a (2AD 5 z l fi m z 2 (m h 2 - m z 2 ) s (5 2 - a 2 ) 

+ D 5 Z % 5 (32m/ { s -5a) + s 2 {5 2 -a 2 f 

- 2m z 2 s(5-a) 2 (2s + 5 2 + 65a + a 2 ) 
+ Am h 4 {s + 5a) (-8m z 2 + {5 + a) 2 ^j 

- Am z 4 {-{5a {5 2 + 18 5 a + a 2 )) + s {7 5 2 + 65 a + 7 a 2 )) 
+ 2m h 2 (?>2m z 4 5a + s{5 + a) 2 (-2 s + [8 + a) 2 ^ 

+ 4m z 2 {2s 2 -5a {5 2 + 10 5 a + a 2 ) + s (3 5 2 + 2 5 a + 3 a 2 ))))))) ,(3.159) 

<f (2) = I (^M {Sm z 2 s + As 2 -s5 2 + sa 2 -A5 2 a 2 ) 

+ c zh,v ( 8m z 2 s - As 2 - s 5 2 + s a 2 + A5 2 a 2 )^j , (3.160) 

G z'h (1) = ^ (-4 Czm m z 2 (2 m h 2 -2m z 2 + s) 

+ Cz/,,1 K 2 - m z 2 - s) {5 2 - a 2 ))) , (3.161) 

Gf^ 0) = - J-, (C ZM (<5 2 - a 2 ) (8 mft ^ ff ( s + ^)-8 m / (6 S 2 + ^a- 5 2 a 2 ) 

+ s 2 ( 5 4 + 6(y 2 CT 2 +fJ 4_ 4s ( 5 2 + fJ 2^ 

+ Am z 2 s (As 2 -5{5-a) 2 a + s {35 2 + 25 a + 3a 2 )^j 

- Am, 2 (-( s {8-a) 2 {s + 5a))+Am z 2 {2s 2 + 5 2 a 2 ))) 
+ 4C ZM (l6m/5a (-s + 5 a) + s {5 2 - a 2 ) 2 {-s 2 + 5 2 a 2 ) 

- Am z 4 ( c [5 2 {5-a) 2 a 2 )+2s 2 {5 2 - 35a + a 2 ) 
+ 2s5a {5 2 -5a + a 2 )) 

+ Am z 2 {2s 3 {5 2 -5a + a 2 ) - 5 3 a 3 {5 2 + 6 5 a + a 2 ) 

+ s5 2 a 2 {75 2 + 25a + 7a 2 ) - s 2 {5 4 - 5 3 a + W5 2 a 2 - 5 a 3 + a 4 )) 

- Am h 4 (-Am z 2 s 5 a + s 2 {5 + a) 2 - 5 2 a 2 (Am z 2 + (5 + a) 2 )) 
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+ 4m h 2 ((<5 + a) 2 (s - 5 a) 2 (s + 5a)+8m z * (s 2 - S 2 a 2 ) 

+ m z 2 (-8s 3 -2b 2 a 2 (5 2 + a 2 ) +2s5a (5 2 + 6 a + a 2 ) 

+ s 2 (5 5 2 -ASa + 5a 2 ))))) , (3.162) 

where we have used the following coupling functions 



D S z \i= 5D^ k [ + aD { ^\ (3.163) 

D S z\ 2 = SD^-aD^p, (3.164) 

D S z\ 3 = ^zty + ^ { ZH k L (3-165) 

D S z a h ,4= ^zty-°V { ZH k L (3-166) 

D S z\ 5 = SD^ + SD^, (3.167) 

D S z\e= SD^-SD^, (3.168) 

D Zh j= vnW,v+° D Zh,A, (3-169) 

D Zh , 8 = <rD { j$-° Di z$t> (3-170) 



Czh,2 
Czh,3 
CzhA 



— o 



[ij]kl ^[ijjkl 
Zh,A ^Zh,V ' 
^[ij]kl _ ^[ij]kl 
^Zh,A ^Zh,V ' 



c 



{ij}M _|_ £,{ij}kl 



Zh,A 



Zh,V ' 



^Zh,A ^Zh,V ' 



(3.171) 

(3.172) 
(3.173) 
(3.174) 



n {ij}fc« 



n [ij]kl 
u Zh,A 



D 



{ij}kl 

Zh,V 



D 



[ij]ki 

Zh,V 



x°A z f^W z 



A 



+ 



°-°*\* X°X> { X °X? h 





a 



V 



(3.175) 



(3.176) 



(3.177) 



(3.178) 
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UJ 



(A-Z) 
Zh 



Re 



x 



C 



QZZh(J X i 



2m% 



s — m A + irriA^A ) s — m z + imz^z 



m h A + (m z 2 - s) 2 -2m h 2 (m z 2 + s)) 



Higgs (A)-neutralino cross term: 



UJ 



(A-x°) 



1 4 



8mis , 
z fc=i 



v°v°4 

(jZhA(j^ x j A 
s — m 2 A + imA^A 



^A,Zh 



4s 



-2 (m^ 2 + m^ 2 — s) s + (m^ 2 — m^ 2 — s) <5 2 ) 



+ (s (-8 (2m h 2 m z 2 - m^ 2 (m h 2 + m z 2 - s)^j s 

+ 2(-((m£-m|) (2m x o 2 + m /l 2 -3m z 2 )) 

+ (^2 m x o 2 + m/, 2 + 3 m z 2 ) a) 5 2 - + m z 2 + s) c5 z 



- 2 (2m x o (m h 4 + (m z 2 -s) 2 -2m h 2 (m z 2 + a)) (s - <5 2 ) 
+ (m£ - m|) <5 (2 (m ft 2 + m z 2 - s) s + (-m^ 2 + m z 2 + s) 5 2 )) cr 
+ (2 s (-(m h 2 -m z 2 f + {m h 2 + m z 2 ) s) 
+ (2(m fc 2 -m z 2 ) 2 - (3m h 2 + 5m z 2 ) s + a 2 ) J 2 ) a 2 ) 
x ' 



s,a, 5, m 2 h ,m 2 z ,m x ^ 



G^fj^ = - <5 ( 4 (m/j 2 - m z 2 - s) s [2m x o - a) + ( 8m x o 3 (m^ 2 - m z 2 
+ 4m x o 2 s (-mh 2 + m z 2 + s) a 

— 2 m x o (m^ 2 — mz 2 — s) (2 mz 2 (s — 5 a) + 2 (s + 5 a) 
+ s (-2s + 5 2 +a 2 ))+a (2m h 4 5 (5 + a) + s 2 (5 - a) (5 + a) 
+ m z 2 s (As - 55 2 + 25 a - a 2 ) + 2m z 4 (-2s + 5 (5 + a)) 
+ m h 2 (-(s (35 -a) (5 + a)) + 4m z 2 (s - 5 (5 + a))))) 
x T s,a,5,m 2 h ,m 2 z ,m x o 2 ) , 

where we have used the following coupling functions 



(3.185) 

s) s 



(3.186) 



c 



[ij]fc 
Zh 



r X°A Z r X° k X°H xlx^Z X»xlH _ 



(3.187) 
(3.188) 



Z-neutralino cross term: 



(z-x°) 



Zh 



fc=i 



[\s-m 2 z + imzVz) I6m 4 s 2 ^ ^ + z - a 



_i_ /-«+ f~p zhs (~Pzh,S 
^Zh,S^Z,Zh ^ZhS Z,Zh 



(3.189) 



G z Z zC = ~ { S (~ 4s [Am^ 2 m z 2 s + 6m z 4 s - Am z 2 s 2 + 2 s 3 - m z 2 s 5 2 

— 4 m^o (m/j 2 — mz 2 — s) (2 mz 2 — s) a + 2 mz 4 5a — mz 2 s a 2 

— 2m h 2 (s 2 + m z 2 (-s + 5a))) + [~lQm x ^m z 2 s 2 
+ 16m x o 3 (rrih 2 — mz 2 — s) (2m/ — s) sa 

— 4 m o (m^ 2 — mz 2 — s) (2 mz 2 — s) a (2 mz 2 (s — 5 a) + 2 m^ 2 (s + 5 a) 
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+ s (-2 s + 5 2 + a 2 )) -8m x o 2 s (s 2 (-m h 2 + s) + 2m z 4 (s + Sa) 

- m z 2 (s 2 + 2m h 2 5a + s (5 2 + a 2 ))) 

+ s {-2m h 2 (m h 2 -s) s (5 + a) 2 + 8m/ (4s -a (6 + 3a)) 

- m z 2 (s(5-a) 2 (-2s + (<5 + a) 2 ) 

+ Am h 2 [As 2 + 5a(5 + a) 2 -s (5 2 +a 2 ))) 

+ m z 4 (s (-6<5 2 + 4 5a - 30a 2 ) + 4a (2m ft 2 (6 + 3a) 

+ 5 (5 2 + 65a + a 2 ))))) f[s,a,5,m 2 h ,m 2 z ,m x o 2 ^ , (3.190) 

= - 2 s 5 (-4 s (-6 mz 4 + (m h 2 - s) (s - a 2 ) + m z 2 (s + a 2 ) ) 

- (-7m z A s5 2 + m z 2 s 2 5 2 - 12 m z 6 5 a + 2m z 4 s5a - 2m z 2 s 2 5 a 

+ 12 m z 6 a 2 - 23 m z 4 s a 2 + 5 m z 2 s 2 a 2 + 26 m z 4 (5 2 a 2 - 5 m z 2 s 5 2 a 2 
+ s 2 5 2 a 2 + 2m z 4 5 a 3 + 2m z 2 s 5 a 3 - m z 2 s a 4 - s 2 a 4 

- m h 4 (5 + a) (-25a 2 + s (5 + a)) 

- 4m x o 2 s (Qm z 4 — (mh 2 — s) (s — a 2 ) — m/ (s + a 2 )) 

+ m fe 2 (l2m z 4 (2s + (<5-a)a) + s(<5 + a)(a 2 (-3<5 + a) + s(<5 + a)) 
+ m z 2 (-8s 2 -45a 2 (5 + a) + 2s (5 2 + 3 a 2 )))) 

x J 7 s,(j,iJ,m^m|,m x o 2 ^ , (3.191) 

G^zvf = 4 s (8 m x o 3 m z 2 s + 4 m x o 2 m z 2 s a 

- 2 m x o mz 2 (2 mz 2 (s — 5 a) + 2 mh 2 (s + 5 a) 

+ s (-2s + 5 2 + a 2 )) - a (-2s 3 - 2m z 4 (3 s + 5 a) + m z 2 s (As + 5 2 +a 2 ) 
+ 2m h 2 (s 2 + m z 2 (-s + 5a)))) + (^32 m x o 5 m z 2 s 2 + 16m x o 4 m z 2 s 2 a 

- 16m x o 3 m z 2 s (2m z 2 (s - 5a) + 2m h 2 (s + 5 a) + s (-2s + (5 2 + a 2 )) 
+ 8m x o 2 sa (s 2 (-m h 2 + s) +2m z 4 (s + 5 a) 

- m z 2 (s 2 + 2m h 2 5a + s (5 2 + a 2 ))) 

+ s a (2 m h 2 (m h 2 - s) s (5 + a) 2 + 8m/ (-4 s + 5 (3 5 + a)) 

- 2m z 4 (4m h 2 5 (35 + a) + s (-15 5 2 + 25 a - 3a 2 ) 

+ 25a (5 2 + 65a + a 2 )) + m z 2 [s(5-a) 2 (-2s + (5 + a) 2 ) 

+ Am h 2 (As 2 + 5a(5 + af -s (<5 2 + <r 2 )))) 

+ 2m x o (8m z 6 (5(j (-s + 5 a) + m z 2 s 2 (5 4 + U5 2 a 2 + a 4 - 8 s (5 2 + a 2 )) 

+ 2 s 3 (-2<5 2 a 2 + s (5 2 + a 2 )) 

+ 2m A 4 (4m z 2 (5 2 (j 2 - 2s5a (-2m z 2 + 5a) + s 2 (c5 2 + a 2 )) 

+ 2m z 4 s (-16s 2 -25a (5 2 - 5 a + a 2 ) +s (7 5 2 + 45 a + 7 a 2 )) 

+ 4m/, 2 (m z 2 s<5a (-2 s + 5 2 + cr 2 ) 
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+ 4m z 4 ( s 2 -5 2 a 2 )- S 2 (-25 2 a 2 + s (5 2 + a 2 ))))) 
x T 



^Z.Zh 



s,a, 5, m 2 h ,m 2 z ,m x o 2 



2s (As (-6m/ + (mh 2 



(3.192) 



s s 



4m x o 2 s (6 m/ 



s) [s 

2 „ , „2 



5 2 ) +m z 2 (s + <5 2 )) a 
i 2 )-m z 2 (s + <5 2 )) a 



m z \s 

2 2 



+ 2m x os (m h 4 + 13 mz 4 - 2 m z 2 s + s 2 - 2 m h 2 (m z 2 + s)) (<T - cr 2 ) 

+ a (12 m z 6 5 (5 - a) + s 2 5 2 (-5 2 + a 2 ) 

- m h 4 (5 + a) (-25 2 a + s (5 + a)) 

+ m z 4 (2<5 2 <r (5 + \3a) + s (-23(5 2 + 2 5 cr - 7a 2 )) 

+ m z 2 s (s (55 2 - 25a + a 2 ) - 5 2 (5 2 - 2 5 a + 5 a 2 )) 

+ m h 2 (l2m z 4 (2s + 5 (-5 + a)) + s (5 + a) (5 2 (6 -3a) + s (5 + a)) 



+ m z 
x T 



-8s 2 - 45 2 a (5 + a) + 2s (35 2 +a 2 ))))) 



X 2 2 2~|\ 

s,a,6,m h ,m z ,m x o j , 
where we have used the following coupling functions 



(3.193) 



^Zh,S 
^Zh,A 



^XiX^Z {ij}k , „XiX°jZ r {ij}k 

a ^zh,A ^ u y ^zuy ' 



c 



XiX°jZ [ij]k I , 



(3.194) 
(3.195) 



yP 

Ai A.j 



c 



c 
c 



{ij}k 
Zh,A 

[ij]k 



Zh,A 
{ij}k 



C 



Zh,V 

[ij]k 
Zh,V 



r X°A Z r X° k X°r xlxlZ X°X° k r 

°A °5 ~r°A °S 



A*l Z C xlx\r 



n x\x ( l z n XjX h T 
A °S 



U T/ -+- L^y L^q 



r X*A Z r xlx»r 



x° fc x^ c x°x° fc 



ZA 



(3.196) 
(3.197) 
(3.198) 
(3.199) 



Contributions to come from s-channel Higgs boson exchange, t- and u-channel neutralino 
exchange and cross terms 



_ r ,(h,H) ~(x°) (mt-x°) 

W X?X?-2A - W ZA + W ZA + W ZA 



(3.200) 



S-channel CP-even Higgs boson (h,H): 



ZA 



=/l,.ff 



^AZr£<X;Xj A 



m?, + im r T r 



(s - a 2 ) (s 2 - 2 (m| + m^) s + ( 



m' 



>> 2 ) 



2m| 



(3.201) 
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T- and U-channel neutralino: 



where 



~,(x°) 



ZA 



2m Z CI 
z k,l 



4 



Xfc Xi 



kl 



(3.202) 



l ki 



ZA,{ij} 



(- 



/,r(o). 



m 



J kl — 



D 



Z ( a + $) T 2 + G J ZA il)j l + G ZA { ° )7 + C ZA,{ij} 



(3.203) 



(a + 6) y 2 



K, 



ZA 
kl 



+ Gf/ ] yi + Gf/ ] y Q ) (s, a, 5, m 2 A ,m 2 z , m^m'o) , (3.204) 

( r K,T(2) r r K,T(l) r r K,T(0) r 

r K,Y(2) v r K,Y(l) v r K,Y(0) v \ / _ r 2 2 2 2 x 

(3.205) 



G'zT = Ys { 2D ZAM -mz 2 -s)5a 

+ D ZA,{ij} ( 2m ^ 2 {s-5a) + 2m A 2 (s + 5 a) 

+ s (-2s + <5 2 + a 2 ))) , (3.206) 

+ D ZA,[ij] ( m ^ 2 ~~ ~~ s ) ^ a ( 2m z 2 (s - 5 a) 

+ 2m A 2 (s + 5a) + s (-2s + 5 2 +a 2 ))) 

+ D+ A{ij} (8 m z 4 (5ct(-s + (5ct) + 8 mV <J <r (s + 5 a) 

+ 477U 2 (s + 5a) (s(5-a) 2 

+ 4m z 2 (s-J^+s 2 (5 i + 68 2 a 2 + a i -As {5 2 + a 2 )) 

- 4m z 2 s (^8 s 2 + 5 (<5 - a) 2 <r - s (3<5 2 + 2<5a + 3<r 2 )))) , (3.207) 

^ ( ° ) = - I ^(24C- Aifo . } ^ 2 S 2 (5 + .) 2 

+ C+ A{ij} (8m z 4 5a (s + Sa) + 8m A 4 5a (s + 5 a) 

+ Am A 2 (s + 5a) (s (5 - o) 2 + Am z 2 (s - 5 cr)) 

+ s 2 (5 4 + 65 2 a 2 + a 4 -4s {5 2 + a 2 )) 

- Am z 2 s (8s 2 + 5{5-a) 2 a-s (3 5 2 + 25a + 3a 2 )))) , (3.208) 
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Gf A W = -y s ((6 + a) (-GD- Mij} m z *s 

+ D+ A {ij} (2m A 2 5a-2 m z 2 (s + 5 a) + s (5 2 + <r 2 ) ) ) ) , (3.209) 
G J Z T = ^ ((* + a) (24m z 2 (- (D- A>[y] (m A 2 - m z 2 + s)) 

+ ^Afol + m z 2 + *)) <5<r 

+ 12L> ZAfo . } m z 2 (2<5cr (m z 2 -25a) -2m A 2 (2s + 5a)+s (5 2 + a 2 )) 
+ ^| Afe} (™z 4 (-32 s + 8 <5 a) + (5 + a) 2 (s (6 - a) 2 + 4 m A 2 5 a^j 
+ 4m z 2 (s (5 2 + a 2 ) -5a (2m A 2 + 5 2 + 65 a + a 2 ))))) , (3.210) 
Gf A {1) = ~\ ((3C ZAm m z 2 + C+ Am ( mA 2 + 2m z 2 -s)) (5 + a)) , (3.211) 

Gf/ ] = ii ((^ + ^) ( 24 ( s (-^ 2 + - * + * 2 ) 

+ (s - 2 5 2 ) a 2 ) + C+ A (8 m A 4 (s + 5 a) - 8m z A (2 s + 5 a) 

+ s ^ 8s 2 + ( 5 2_ fJ 2^_ 4s (j2 + ff 2^ 

- 4m z 2 (6s 2 + 5a (<5 2 + 6 <J a + a 2 ) -s (3 5 2 + 25 a + 3a 2 )) 

+ 4m/ (6m z 2 s-(s + 5a) (4s - (<5 + ^) 2 ))))) , (3-212) 

" ^ D zA M sSa + D zA, {i j} ( s ~ U2 ) "*) > ( 3 - 213 ) 
Gjf (1) = (2 , (-6 D^ >W} m Z 2 s(5 + a) 2 

- D+ A[ij] 5a (8m z 4 + m z 2 (6 s - 5 2 - 6 5 a - a 2 ) 

+ 2s (-s + 5 2 + o- 2 ) + m A 2 (-8m z 2 + 2 s + 5 2 + 65 a + a 2 ))) 
+ D+ A{ij} (\6m z A s (s-5a) + 2m A 2 (s 5 (5 - a) 2 a - A5 3 a 3 
+ 8m z 2 s (s + 5a) - s 2 (5 2 + 4 5 a + a 2 )) 

- 2m z 2 (8s 3 -45 3 a 3 + s5a (5 2 -W5a + a 2 ) - s 2 (5 2 + 4 5 a + a 2 )) 
+ s (_ 2s 2 (j 2 + a 2 ) -4<5 2 <r 2 («5 2 + a 2 )+s («5 4 + 65 2 a 2 + <7 4 )))) , 

(3.214) 

Gff (0) = ^ {D + ZAm (64 mz 6 Sa (s - 5 a) 

+ 8m A 4 5a (s + 5 a) {-8m z 2 + (<5 + a) 2 ) 

_ s (5 2 -a 2 f (-A5 2 a 2 + s (5 2 + a 2 )) 

+ 8m z A ( K As 2 (5-a) 2 + s5a (5 2 - 14 5 a + a 2 ) 
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+ 5 2 a 2 (5 2 + W5a + a 2 )) - Am z 2 (s 2 (5 - a) 2 (5 2 + 4 5 a + a 2 ) 

+ A5 3 a 3 (5 2 + 65a + a 2 ) - sSa (5 A + 12 5 3 a + 6 5 2 a 2 + 12 5 a 3 + a 4 )) 

- 4m A 2 (32m z 4 (s 2 - <5 2 a 2 ) - (S + <r) 2 (45 3 a 3 - s 5a {5 + a) 2 

+ s 2 ( ( j2 + <T 2 ))-4mz 2 (s 2 (<5 2 + 4<5a + <7 2 )-s<5<7 (<5 2 + 6 5 a + a 2 ) 

- 5 2 a 2 (5 2 + 65a + a 2 )))) 

+ 25 a (24 (l?-^.., m z 2 (-m A 2 + m z 2 ) s (5 + a) 2 

+ D+ AM (32 m z e (s-5a) + S 2 {5 2 - a 2 ) 2 

- 2m z 2 s(5-a) 2 (2 s + 5 2 + 6 5 a + a 2 ) 
+ 4m A 4 (s + ^) (-8m z 2 + (<5 + (j) 2 ) 

- 4m 2 4 (-(<5<r {5 2 + 185a + a 2 )) + s (7 5 2 + 6 5 a + 7 a 2 )) 
+ 2m A 2 (32m z 4 5a + s(5 + a) 2 (-2 s + {5 + a) 2 ^ 

+ Am z 2 (2s 2 -5a (5 2 + 10 5 a + a 2 ) + s (35 2 + 25a + 3a 2 )))))) (3.215) 

G K ZA (2) = ^ (" {c + ZMij} s{5 + af) +±C ZMi3} (2m z 2 s - s 2 + 5 2 a 2 )) , (3.216) 

G Z 'I W = j- g (Sa (±C- A [ij] m z 2 (2m, 2 -2m z 2 + s) 

+ C %am ~ ™z 2 -s)(5 + a) 2 )) , (3.217) 

Gfr ] = 7i^{-{c + ZAm (5 + a) 2 (8m A Ua(s + 5a) 

- 8m z 4 (6s 2 + s5a-5 2 a 2 ) + s 2 (<5 4 + 6 5 2 a 2 + a 4 - 4s (<5 2 + a 2 )) 
+ Am z 2 s (4s 2 -5{5-a) 2 a + s (35 2 + 25a + 3<r 2 )) 

- 4mA 2 (- (s(5-a) 2 (s + <5a)) +4m z 2 (2 s 2 + <5 2 a 2 )))) 

+ *C- At{ij} (l6m z 6 5a {-s + 5a) + s(5 2 -a 2 ) 2 (- s 2 + 5 2 a 2 ) 

- Am z 4 (-(«5 2 (<5-a)V)+2s 2 (5 2 -35a + a 2 ) 

+ 2s5a (5 2 - 5a + a 2 )) + Am z 2 (2s 3 (5 2 - 5 a + a 2 ) 

_ S 3 a 3 ( 5 2 + 66a + a 2^ +s§ 2 a 2 ( 7S 2 + 2 5a + 7a 2 ) 

_ s 2 (s±-6 3 a + 105 2 a 2 -5a 3 + a 4 )) 

- 4m A 4 ( K -Am z 2 s5a + s 2 {5 + a) 2 -5 2 a 2 (^m z 2 + (5 + af^j 
+ 4m A 2 ((5 + a) 2 (s - 5 a) 2 (s + 5a) + 8m z 4 (s 2 - 5 2 a 2 ) 

+ m ^ (_ 8s 3_ 2(5 2 a 2 ( J 2 +(T 2^ ) 

+ 2s(5fj (<5 2 + (5a + a 2 ) +s 2 (55 2 -45ct + 5c7 2 ))))) , (3.218) 
where we have used the following coupling functions 
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ZA 



I 4 2 

8m 2 7 s ^ ^-f Re 

Z l=i r =l 



(jZrA(j x< i x< j r 

s — m 2 + im r T r 



, (3.235) 



U r,ZA 



^r,ZA 



5 (— m^ 2 + m^ 2 + s) cr + |4m x o 2 s (— niA 2 + m z 2 + s) a 

- 4m x o (m A A + (m z 2 - sf -2m A 2 (m z 2 + s)) (s - cr 2 ) 
+ a (2 m/ cr (5 + cr) + m z 2 s (4 s - 5 2 + 2 «J a - 5 a 2 ) 

+ s 2 (-<5 2 +cr 2 ) +2m z 4 (-2s + a (5 + cr)) 

+ m A 2 (s (5 2 - 25a -3a 2 ) + Am z 2 {s - a {5 + a))))) 

x T s^a^^m^^m^^m^ 2 ^ , (3.236) 
-4s [m A 2 (2 s + (2 m x o - a^j a^j + s (-2 s - 2 m x o cr + a 2 ^ 

+ m z 2 ^2s-2m x ocj + cj 2 ^ - (-2m z A s 5 2 + 2m z 2 s 2 5 2 
+ 8 ?^ x o 3 {mA 2 — m z 2 — s) s a + 4 m^ 4 s 5 a — 4 m^ 2 s 2 5 a 

- 6m z 4 s a 2 + 6m z 2 s 2 a 2 + 2m z 4 5 2 a 2 - 5m z 2 s 5 2 a 2 + s 2 5 2 a 2 
+ 2m z 4 (5a 3 + 2m z 2 S(5CT 3 -m z 2 scj 4 -s 2 o- 4 

- 2m A 4 (5 + a) {- {5 a 2 ) + s (5 + a)) 

+ 4m x o 2 s (m A 2 (2 s -a 2 ) + s (-2s + a 2 ) +m z 2 {2s + a 2 )) 

- 2 m x o [mA 2 — m z 2 — s) a (2 m z 2 (s — 5 a) + 2 mA 2 (s + 5 a) 

+ s (-2s + 5 2 + a 2 ))+m A 2 [s (5 + a) {a 2 {-35 + a) + 2s (5 + a)) 

- Am z 2 (4s 2 + 5a 2 {5 + a) - s (5 2 + 2a 2 )))) T s,a,5,m 2 h ,m z ,m x o 2 , 



(3.237) 



where we have used the coupling functions 



Cf = 



(c^ A y + (ct- z y (c. 



(3.238) 
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A 



(3.239) 



x?x? -> WW 



Contributions to cl> come from s-channel Z and Higgs boson exchanges, t- and u-channel 
chargino exchange and cross terms 



w v o_o 



(h,H) , _(Z) , -(**) , , M-X 



+ to- 



WW 



+ LO 



WW 



(3.240) 



S-channel CP-even Higgs boson (h,H): 



~(h,H) 
t^WW 



E 

r=h,H 



(jWWrfjXiXjr 



s — m 2 + im r T r 



(s-a 2 ) (s 2 -Asm w + 12m^) 



8m^ 



(3.241) 



S-channel Z-boson: 



LO 



WW 



C 



WWZQXiX] 



4 



s — m z + imz^z 



(s - a 2 ) (2s + 5 2 ) (s 3 + 16s 2 m 2 v - 68s - 48m^) 



24s 



(3.242) 



T- and U-channel chargino: 



lo 



(x±) _ 1 



WW 



2 



m+m + 



x 



-+- to + -+- ic^ 

xt kl kl 



(3.243) 



where 



J fci — 



J kl — 



,y 

= (3D Sa 



w \v'To + C WWD2 [ s 



(s - Am 2 w ) y 2 



+ Gww^o) ( s i a i^ m wi m W' m l±' m l±) ' 
tw^w % + G^'Ti + G J ™% 

+ D WW,C,2 (* " "4) 3>2 + 



-\ i \ / r 2 2 2 2 ■ 

•i, , K 3^0 {s,cr,5,m w ,m w ,m +,m +, , 

/ A-fc A/ 

(-Cw,d^4 - C++,,d (s - a 2 - «5 2 - 2m^) T 3 



\ .... . . . . 

r K,T(2) r r K,T(l) r r K,T(0) r r+ - y 
+ U WW 1 '2 + { - r WW J 1 + L 'WW J + ^WW,D,2 ^ 

+ G ww 2) y^ + ^W^o) (a, <r, m^, m 2 w , m 2 ±, m 2 ±) , 



(3.244) 



(3.245) 



(3.246) 
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G ww = \ C ww,d (16 m^ 4 - s (5 2 + a 2 ) + Am w 2 {-2 s + 5 2 + a 2 )) , (3.247) 

G ww = ( C ww,d (-64 mw 6 + 64 m w 4 s - 40 m w 4 5 2 + 16 m w 2 s 8 2 

- 4 m w 2 5 4 + s S 4 - 40 m w 4 a 2 + 16 m^ 2 s cr 2 - 8 m^ 2 <5 2 cr 2 

- 2s5 2 a 2 - Am w 2 a 4 + s a 4 ) 

+ Cw,D (72 m^ 4 (S 2 - a 2 ))) , (3.248) 

G 1 /^ = Y% (C ww ^ 2 (-64 m w « + s(5 2 - a 2 ) 2 + 8m w 4 (5 2 + a 2 ) 
+ Am w 2 (2s-5 2 -a 2 ) (5 2 + a 2 ) 

+ 4C W,D,2 (6 m w 4 + 4 m w 2 s - s 2 ) (5 2 - a 2 ) ) ) , (3.249) 

C = - \ m lv (d { w + } 5 (5 2 - 2m 2 w ) + C { w ^a {a 2 - 2m 2 w )) , (3.250) 

Gw7 ] = H2 ™ w 4 S <5 + 4 m w 2 S <5 3 + s <5 5 

+ 12mvi/ 2 S(5CT 2 - 8m w 2 S 3 a 2 - s5a 4 ) 

+ (-32 mvi/ 4 s cr + 12 m^ 2 s5 2 a-s5 4 a + A m w 2 s a 3 

- 8m w 2 5 2 a 3 + sa 5 )) , (3.251) 

G W) = Ok (D^ 2 (l2m 2 w -8s + 3(5 2 + a 2 )) 

+ 3 D ww,c,2(° 2 -6 2 )) , (3-252) 

G ^ 0) = ib (-^,2 ( 16 ^ 6 s - * 2 (* 2 - - 2 ) 2 

- 8m^ 4 (lis 2 + 5 2 cr 2 - 5s (<5 2 + a 2 )) 

+ m w 2 s (lQ s 2 + 5 4 -U5 2 a 2 + a 4 -2s {5 2 + a 2 ))) 

+ 6sm 2 w D% WC2 (4m w 2 - s) {5 2 - a 2 )) , (3.253) 



G^ 2) = — (C++, C (-40 m w 4 s + 32 m w 2 s 2 - 12 m w 2 s 5 2 + 4 s 2 <5 2 - 3 s 5 4 

- 12 m w 2 sa 2 + As 2 a 2 + 8 m^ 2 <5 2 a 2 - 2 s 5 2 a 2 - 3 s a 4 ) 

+ Sm^tfV) , (3.254) 

G ™ = ( 128 mw " s ~ s {s ~ 52 ~ a2) {52 ~ a2) 2 

- 16 m w 4 (4 s 2 - 5 2 a 2 ) - 2 m w 2 (8 s 2 (5 2 + a 2 ) + 4 5 2 a 2 (5 2 + a 2 ) 

- s(35 4 + 25 2 a 2 + 3a 4 ))) + 72m^s (,5 2 - a 2 ) C+~ WD ) , (3.255) 

G ™ = " 2W ( 1024 mw% s2 ~ 256 * 2 ^ " 48 mH/4 S * ^ 
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+ 8 m w 2 s 2 5 6 + s 2 5 8 - 256 m w 6 s 2 a 2 

+ 512 m w 6 s5 2 a 2 - 416 m w 4 s 2 5 2 a 2 - 64 m w 4 s S 4 a 2 - 8 m w 2 s 2 5 4 a 2 

- 16m w 2 sS 6 a 2 -4s 2 5 6 a 2 -A8m w 4 s 2 a 4 -64m w 4 s5 2 a 4 

- 8 m w 2 s 2 5 2 a 4 + 128 m w 4 S 4 a 4 
+ 32m w 2 s5 4 a 4 + 6s 2 5 4 a 4 

+ 8 m w 2 s 2 a 6 - 16 m w 2 s5 2 a 6 -4 s 2 5 2 a 6 + s 2 a 8 ) ) , (3.256) 

G ww 2) = (tw,D,2 ( 8 m ^ s " 32 mw<2 * 2 + s 64 + 8 mw/2 ^ ^ " 2 s 52 a2 + s ff4 ) 

+ G W,D,2 ( 8 s ^ - 2 s 2 ^ - 8 m w 2 s a 2 + 2 s 2 a 2 ) ) , (3.257) 

G ^ 0) = 2k (" 4C ™ 2 ^ ^-^s + s 2 (J 2 -a 2 ) 2 

- 8mvK 4 (6s 2 + 105 2 a 2 -s (5 2 + a 2 )) 

+ 4m w 2 s (-J 4 + 6<5 2 a 2 -a 4 + 2s (5 2 + a 2 ))) + C+~ WjDj2 (s (<5 2 - <x 2 ) 4 

+ 16m w 2 (<5 2 -a 2 ) 2 (-4 s 2 + 5 2 a 2 ) 

- 256 mvK 6 (8 s 2 - 5 2 a 2 - 2 s (<5 2 + a 2 ) ) + 16 m^ 4 (48 s 2 (<5 2 + ct 2 ) 

+ 20 5 2 * 2 (5 2 + a 2 )-s (9 5 4 + ll8 5 2 * 2 + 9a 4 )))) , (3.258) 



where we have used the following coupling functions: 



D' 




(3.259) 
(3.260) 
(3.261) 
(3.262) 



WW,D,2 — 



WW,D — 

<+± _ 
WW,D — 




(3.263) 
(3.264) 
(3.265) 
(3.266) 




(3.267) 
(3.268) 
(3.269) 
(3.270) 
(3.271) 



<{-+} _ r<-+ j- r<+- 

WW,2 — ^WW,2 ^WW,2 ' 
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a 

D 
D 



+} 

WW,2 

{-+} 
WW,2 

+} 
WW,2 



r<—+ _ fi-\ — 
U W W,2 ^WW,2 ' 



D 



WW,2 + D 



WW,2 ' 



D~ + - D + ~ 
^WW,2 ^WW,2 ' 



(3.272) 
(3.273) 
(3.274) 



D ± 
^ww 

^WW,2 



D 



± 

WW,2 



Re (c wwk {pww,l 
Re {^WW,k (Pww,l) ) ' 

Re (pww,k^ww,i) ' 
Re {pww^ww^i) > 



(3.275) 
(3.276) 
(3.277) 
(3.278) 



a 



WW 



^ww 



D 



WW 



^WW,2 
^WW,2 



D 



WW,2 



D 



Re (C WWk [C wwl 
Re (pww,k {pww,l 



D WW — Re { D WW,k ( D \VW,l) ) ' 



WW,2 



Re I R>ww,k Pwiu 



Re {pWW,kPwW,l) ' 
Re {pWW,kPwW,l) ' 
Re {pww,k®wwj) ' 
Re {^WW,k^WW,l) ' 



(3.279) 
(3.280) 
(3.281) 
(3.282) 
(3.283) 
(3.284) 
(3.285) 
(3.286) 



^WW,k — Ul 



xix jW- 

V 



a 



xtx°w- 



V 



n ± _ r xtx°w- ( x +x?w- 
u ww,i ~ a I °y 



± c xt^w- 



±c xtx°w- 



C 



C 



xtx?w- 

A 

xtx'lw- 



(3.287) 
(3.288) 



Higgs (h, H)-chargino cross term: 



{^H-xt 



WW 



5> 



k=l 
C: 



TT s — m$ + im r r r 

,r=h,H 



E 



x I ^WW^r.WW ^WW,k U r,WW 



(3.289) 
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(~PwW,k 
U rWW 



/^WW,k 

^r,WW 



Am +s (a (2m w 2 -s) s+ Um + 2 (2m w 2 - s) s + 8m w 4 (2s -3a 2 ) 

- 2m w 2 s (5 2 + 25a -3a 2 ) + s 2 (5 2 + 26a - a 2 )) 

x T s.a^.m^.m^.m^ ) , (3.290) 
a (As (l6m w 4 + 4m x + 2 (2m w 2 + s) +s (2s- (5 + a) 2 ^ 

- 2m w 2 (2s + (5 + a) 2 ^ + (-Um w & s + lQm x + 4 s (2m w 2 + s) 



+ 8m w 4 (As 2 + s(S-a) 2 -2b 2 a 2 

+ s 2 (5 4 -8s5a + A5 3 a + 25 2 a 2 + A5a 3 + a 4 ) 

+ 2m w 2 s (5 4 + A5 3 a + 1A5 2 a 2 + A5a 3 + a 4 

+ s (-65 2 + A5a - 6a 2 )) + 8m x + 2 s (Am w 4 + s (2s-(5 + a) 2 ^ 



- 2mw 2 (s + (6 + °") 2 ))) F s, a, 6,171^ ,m^/,m 
where we have used the following coupling functions 



(3.291) 



r ± r xtx { ]w- ( x +x°w- 



(3.292) 



Z-chargino cross term: 



WW 



E 

fc=i 



v o v o 7 



+ imzTz ) 96771^5 



J3Z 



x 1 u WW,k Kj Z?WW ^WW,k u Z,WW 



(3.293) 



U Z,WW 



24 m x + s (5 (48 m w 4 - A s 2 + (-48 m w 6 + 2 (48 m w 4 - 4 s 2 ) 



- Am w 4 (lAs + 35 2 + 65a - 17 a 2 ) + Am w 2 s (4s-3a 2 ) 
+ s 2 (5 2 + 25a-a 2 )) T s,o-,5,m^,m^,m x + 2 ) , 

- 4 (-192r7H F 6 s + 48m x + 4 (2m^ 2 -s) s 

+ 6m w 2 s (2As 2 + S 4 + A5 3 a - Q5 2 a 2 + A5 a 3 + a 4 

+ 4s (25 2 -25a -3a 2 )) - s 2 (-8s 2 + 35 4 + 125 3 a + W5 2 a 2 

+ 125cr 3 + 3cj 4 + 2s (5 2 - 65a + a 2 )) 

- 8m w 4 (28 s 2 + 10 5 2 a 2 + s (-7 S 2 + 6 5 a + 11 a 2 )) 



(3.294) 



- 41 - 



+ 24m x + 2 s {Am w 4 + m w 2 (a s - 2 [5 + <r) 2 ) +s (-s + (8 + ^) 2 ))) 

- 3 (256 mp/s + 64m x + 6 (2m w 2 - s) s 

+ s 2 (5 + a) 2 (5 4 - 8 s 5 a + 4<5 3 cr + 2<*V + 45 a 3 + a 4 ) 
+ 32m^ 6 (16s 2 + 25 2 a 2 + s(5 2 + 65a- 5a 2 )) 

- 2m w 2 s (qAs 2 5 a + (5 2 - a 2 ) 2 (5 2 + Q5a + a 2 ) 

+ 4s (2 5 4 - 5 2 a 2 - 10 5 a 3 - 3 a 4 ))- 16 m w 4 (-(5 2 a 2 (5 + a) 2 ^J 
+ 5s 2 (5 2 - 25a + 3a 2 ) +s (5 4 + 2 5 3 a - 11 5 2 a 2 - A5 a 3 - 2a 4 )) 
+ 16m x + 4 s (2ra^ 2 (8 s - 3 ((5 + a) 2 ) + s (-4s + 3 (5 + a) 2 )) 

- Am x + 2 (96m w 6 s + 16m w 4 (5 s 2 - s 5 2 + 2 s a 2 + 5 2 a 2 ) 

- 2m w 2 s (32s 2 + 3(5 4 + 12(5 3 cj + 25 2 a 2 + 12 ( 5CT 3 + 3cj 4 

- Asa (85 + 5a)) + s 2 (35 4 + 12 5 3 a + 1A5 2 a 2 + 125 a 3 + 3a 4 

- As (5 2 +A5a + a 2 )))) T\s,a,5,m 2 w ,m 2 w ,mA , (3.295) 



where we have used the following coupling functions 



D 



± 

WW,k 



C 



V 



V 



(3.296) 



Y P Y 9 

Ai A.j 



zz 



Contributions to come from s-channel Higgs boson exchange, t- and u-channel neutralino 
exchange and cross terms 



^x\x]^ zz ~ VJzz 



LO 



(h,H) , . (X°) , ~(h,H- X l) _ 



+ u zz + u zz 



S-channel CP-even Higgs boson (h,H): 



(3.297) 



to 



(h,H) 
ZZ 



E 

r=h,H 



(jZZrfjXiXj'r 



s — m 2 . + im r T r 



(s - a 2 ) (s 2 -Asm 2 z + 12m 4 z ) 



16m z 



(3.298) 



T- and U-channel neutralino: 



u zz - 



1 4 

E [ m xl m xfkl Z +rn x ojZ z + 
z k,l=i 



(3.299) 



where 
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+ C z'z,p ( s - 4*4) ^ + G5J (0) ^o) (a, <r, 5, m|, m|, m 2 o,m 2 o) , (3.300) 

Jki Z = (3-1 (C^a + D+ z 5) T 2 + Gg (1) T! + G J ™% 

~ iPtzO + D%~ z 5) (a - m|) y 2 + Gf^y, 

+ G z ' z ^y Q ^ (s,a,5,m%,m z ,m^o,m^o) , (3.301) 

= (-^ P T 4 - (s - a 2 - 5 2 - 2m|) T 3 

r K,T(2) T r K,T(l) r r K,T(0) T 

+ Cz^My^ + Gzz i2) y2 + G z ' z i0) y ) (s,a,5,m 2 z ,m z ,m 2 xl ,m 2 x o), 

(3.302) 

G g {1) = ^(C^p (l6m/-a (5 2 + a 2 )+4m/ (-2 a + 6 2 + a 2 ))) , (3.303) 

G*™ = 1 (c^> (-64m/ + a (<5 2 - a 2 ) 2 + 4m/ (4 a - <5 2 - a 2 ) (S 2 + a 2 ) 

+ 8m z 4 (8s - 5 (5 2 + a 2 ))) + 72 m/ (5 2 - a 2 ) C^T M ) , (3.304) 

4'z (0) = ^ ("64m/ + a (6 2 - a 2 ) 2 + 8m/ (<5 2 + a 2 ) 

+ 4m/ (2s-5 2 -a 2 ) (<5 2 +<r 2 )) 

+ 4 (6 m/ + 4 m/ a - s 2 ) (5 2 - a 2 ) CfjT M ) , (3.305) 

= " \rn 2 z {D+ z S + C+ z a) (-2 m z 2 + a 2 ) , (3.306) 
G J>n°) = M (£)+- 5 (-32 m/ a + 4 m z 2 (a <5 2 + 7 a a 2 - 4 <5 2 a 2 ) 
+ s (-5 2 + a 2 ) 2 ^j +C+ Z a (-32 m z A s 

+ 4m/ (7a<5 2 + aa 2 -4<5 2 a 2 )+a (-<5 2 + a 2 ) 2 )) , (3.307) 
= ~\ (3m| - 2a) (£+/ + C+ z a) (3.308) 

GJ? 0) = ^(^ z -a(l6m/a-a 2 (^-a 2 ) 2 

- 8m/ (a (lla-95 2 ) -2 (a - <5 2 ) a 2 ) 

+ m/a (l6a 2 +5 4 -185 2 a 2 + (j 4 -4a (5 - a) (5 + a))) 

+ D+ z <5 (l6m/a-a 2 (<5 2 - a 2 ) 2 

- 8m/ (lla 2 -2a < 5 2 -9ao- 2 + 2 ( 5 2 cr 2 ) 

+ m/a (l6a 2 +5 4 -185 2 / + <r 4 + 4a (5 - a) (6 + a)))) , (3.309) 
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Gf™ = - \-Czz,P {40mz 4 s + s (35 4 + 25 2 a 2 + 3a 4 -4s (5 2 + a 2 )) 

- Am z 2 (8s 2 + 45 2 a 2 -3s (<5 2 + a 2 ))) , (3.310) 
<f (1) = I ^(c^(l28^ 6 s-s(s-^ 2 -. 2 ) (5 2 -. 2 ) 2 

- 32 m z 4 (2 s 2 - 5 2 a 2 ) - 2 m z 2 (8 s 2 (<5 2 + a 2 ) + 8 <5 2 cr 2 (5 2 + a 2 ) 

- s (3 5 4 + 10 5 2 a 2 + 3 a 4 ))) +72 m z 4 s (5 2 - a 2 ) C^ M ^j , (3.311) 

riK,T(0) 1 ^D,+ / 1n0/) 8 2, 2 / r2 2\4 

zz 256? zi? - p ( 1024m ^ s + s ( 5 " G ) 

- 256m z 6 s (-45 2 a 2 + s (5 2 + a 2 )) 

+ 8m z 2 s(<5 2 -<7 2 ) 2 (-4,5 2 a 2 + s (<5 2 + a 2 )) 

- 16 m z 4 (_l6,5 4 a 4 + 8s<5 2 a 2 (,5 2 + a 2 ) 

+ s 2 (3 5 4 + 26<5 2 <j 2 + 3a 4 ))) , (3.312) 

G K zz m = Ys ( C ^ p ( 8mz4 s + s ( s2 -' j2 ) 2 + m z (" 32s2 + 165V ) 

+ 2C%+ M s (4m| - s) (<5 2 - a 2 ))) , (3.313) 

G ^ (0)= ^(^( s2 (^-- 2 ) 4 -32-z 2 S (5 2 -a 2 ) 2 (2s 2 -* 2 a 2 ) 

- 512 m z 6 s (4 s 2 - 5 2 a 2 - s (5 2 + a 2 )) 

+ 16m z 4 (16<5 4 <t 4 + 48s 3 (5 2 + a 2 ) - s 2 (9 5 4 + 94 5 2 a 2 + 9 a 4 ))) 
+ Czz,M (" 4s2 - ° 2 ) (32m z 6 + s (J 2 - <r 2 ) 2 

+ Amz 2 {2s-S 2 -a 2 ) (5 2 + a 2 ) + 8m z 4 {-6 s + S 2 + a 2 )))) , (3.314) 
where we have used the following coupling functions 



Czz,p — C ZZk \C zzl J +D ZZk [D zzl J , (3.315) 

C zz ± m = C zz,k{ C zz,i) ~ D zz,k{ D zz,i) ' (3.316) 

c zz = c zz,k ( C zz,i) + C zz,k (pzz,i) ' (3.317) 

c zz = c zz,k ( C zz,i) ~ c zz,k (pzz,i) ' (3.318) 

D+ zz = D zz,k ( D zz,i)* + D zz,k ( D zz,i)* > (3-319) 

D zz = D + zz ^ k (D zz i y - D ZZk (D+ Zji y , (3.320) 

C± Ztk = C$* Z (C^ Z ) * + ct° jZ (cf*° Z ) * , (3.321) 
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C% z , = C^ Z (C^ Z ) * ± ct° iZ (C^ Z ) * , (3.322) 
Dl z , k = ct° Z {cf^ Z ) * ± cf^ Z (C^ Z ) * , (3.323) 
flf z>l = cf^ Z (C^ z y ± C^ Z (C^ z y ■ (3.324) 



Higgs (h, H)-neutralino cross term: 



h,H-xl) 



zz 



E E* 

r=h,H k=l 



rtZZrr< 



v°v°r 

Aj Aj 



s — ml + im r T r 



32m z s 



, ^y- f~Pzz,k , /-<+ (~p zz h 
X 1 ^ZZ.fc^r-.ZZ Ly ZZ,fc U r,ZZ 



(3.325) 



y-Pzz,k 
U r,ZZ 



G 



~^ZZ,k 

r,ZZ 



4m x os (4 s (-2m/|s) + (4m x o 2 s (-2m z 2 + s) -8m z 4 (2s -3a 2 ) 
+ 2m/s (<5 2 - 3a 2 ) + s 2 (-5 2 + a 2 )) T s,a,5,m 2 z ,m 2 z ,m x o 2 ) , (3.326) 

- a (4s |l6m z 4 + 4m x o 2 (2m z 2 + s) 

+ s (2s -5 2 -a 2 ) - 2m z 2 (2s + S 2 + a 2 )) 

+ (-64m z 6 s + 16m x o 4 s (2m z 2 + s) + s 2 (5 2 - a 2 f 

+ 8m x o 2 s (4m z 4 + s (2s -5 2 -a 2 ) - 2m z 2 (s + 5 2 + a 2 )) 

+ 8m/ (As 2 -A5 2 a 2 + s (35 2 + a 2 )) 

+ 2m/s (5 4 + 6<5 2 <t 2 + <t 4 -2s (5 2 +3a 2 ))) 



x .F 



s,cr, (5, m|,m|,m o 



where we have used the following coupling function 



^zzm — u y 



yP 

Ai A.j 



ff 



(3.327) 



(3.328) 



Contributions to Co come from s-channel Z and Higgs boson (both CP-even and CP-odd) 
exchanges, t- and u-channel sfermion exchange and cross terms 



,(h,H) , (A) (Z) ~(f) ,M,H-f) 
J ff +UJ ff +UJ ff +U ff 



UJ XiX°^ff = UJ ff +u))f "'+co^ "' +io)f " / +oJ K f f : (3.329) 



.(A-Z) JA-f) (Z-f) _ 
J ff +UJ ff +UJ ff 
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S-channel CP-even Higgs boson (h,H): 



~ (h,H) _ y- 
U ff ~ 

r=h,H 



°5 °5 



s — m 2 + im r T r 



(s - a 2 ) (s - 4m 2 ) ; 



(3.330) 



S-channel CP-odd Higgs boson A: 



U) 



(A) 
ff 



Op O p 



s — m 2 A + imA^A 



s(s-5 2 ) ■ 



(3.331) 



S-channel Z boson: 



(Z) 
ff 



s — m z + imz^z 



3m z 



x 2 



a 



XiXjZ 



+ 2 
+ 
- 6 
+ 
+ 



C 



v 



XiXjZ 



4 

m z 



r ffz 



r ffz 



(s + 2m)) (2s + a 2 ) 



C 



ffz 



(s - 4m 2 ) + 
%(s + 2m}) (s - a 2 ) 
(sm| (s-a 2 ) +m} (3sV - 6s m|a 2 - m% (As- 7a 2 )))^j 



C 



C 



ffz 



XiXjZ 



ffz 

V 



' m%(s + 2m 2 f ) (5a - 25) 



(s m% (5a - 25) + 2m) (m% (A5 - 13a) + 6s m|a - 3s 2 a)) 



a 



XiXjZ 
V 



r ffz 



' m%(s + 2m}) (5a + 46) + 



C 



ffz 



(sm% (5a + 45) 



- 2mj ((3s 2 - 6s m 2 z + 7m z ) (a + 25) + 6m z a))))) ; 
T- and U-channel sfermion: 



(3.332) 



"ff 



in 

a,b 



-{if} 



+ 



1 



( r 2 2 2 

{s,a,d,mf,mj:,m; ,m 



(3.333) 
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Gfp = - 8 m f (2 C% } m f (S 2 - a 2 ) + 2 ro/ (* 2 - a 2 ) 

+ C% j} (l6m/+(5 2 -, 2 ) 2 + 8 m/ 2 (5 2 +a 2 )) 

+ £>? {ij} (l6m/ + (5 2 - a 2 ) 2 + 8m/ (<5 2 + a 2 )) , (3.334) 

<$°> = ^,(W / + (. 2 - 5 2 ) 2 ) 

+ 4 (2A% } m) (2s - a 2 - 5 2 ) + A% } (s - 2m 2 ) (a 2 - S 2 ) 

+ 2m, (s - 6 2 ) - S^d (s - a 2 ) + B^m, (a 2 - <5 2 ))) ; 

(3.335) 

Higgs (h, H)-sfermion cross term: 



CO 



(h,H-f) 
ff 



-YRe 



E 



c{ fr c 



r s-m-t + im T i. r 

r=h,H 



Y [2m f oC» +{ij} [4+ [Am 



Am 4 



+ 2 s — 5 2 — a 2 ) T ^s, a, 5, m 2 , m 2 , mj 2 ^ 
+ D b +{IJ} (4, + (4m ; 2 (s - 2a 2 ) + s (Am f ~ 2 - 5 2 + a 2 )) 
s,a,S,m 2 ,m 2 f ,mf 2 



x T 



(3.336) 



Higgs (A)-Z cross term: 



~(A-Z) u 
U ff =Re 



r ffA r X iX°j A 
^ p p 



U A U A 



x?x° 3 z 



s — m A + imA^A ) \ s — m z + imzTz 



4am 



mi 



l(s-5 2 ) (m 2 z -s)- 
(3.337) 



Higgs (A)-sfermion cross term: 



CO 



-s5> 



r ffA r X° l X° J A 

P P 

-m\ + imA^A 



x {p+{ij}(^ m f ( s ~ ft 2 ) a F s,a,S,m 2 ,m 2 ,mj: a 2 ^ 

+ D lm ( s (" 4 + (" 4 m h + 4 ™/ 2 " 52 + ^ 2 ) 



x T 



X 2 2 2 

s,a,d,mjr,mf,mj a 



(3.338) 
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Z-sfermion cross term: 



(*-/) 



ff 



Am\ 

6 a 



s — m 2 z + im z T z 



C ffZ C ^ Z Z (+)ff + C {r fZ C^ Z Zj-Vf 



(3.339) 



where 



Z Wff 



- AD a +{ij} m f (-3m z 2 + s) a + C^ {ij} m z 2 (-4m /a 2 + 4m / 2 + 2s + i 2 + a 2 ) 

- i (aD^^ m,fa (Arrij (3m z 2 - s) - 4m| s + 4m^ 2 (-3m| + s) 
+ <5 2 ml + 5 2 s + 3m 2 z a — a 2 s) 

- C£ {ij} (16m 4 - m| - 8m 2 , (4m ^ 2 m| + 2ct 2 (s - <5 2 ) + m| (2s + <5 2 - 5ct 2 )) 
+ m 2 z (l6m /a 4 + (5 2 - a 2 ) (4s - <5 2 - a 2 ) - 8m /a 2 (5 2 + c 



x T 



X 2 2 2 

s, a, o,m f ,m f ,m ? 

J •> Ja 



Z (-W = AD a _ M m f 5 + m z 2 (-4m /a 2 + 4m/ + 2 s + 5 2 + a 2 ) 



- - [C a _ {ij} m z 2 (^16m /a 4 + 16m / 4 + 16m / 2 s + 4s5 2 + ( 5 4 -4scj 2 
- 2 5 2 a 2 + a 4 + 8m/ 2 (<5 2 - a 2 ) -8m /a 2 (4m/ + 5 2 + cr 2 )) 
+ AD a _ {ij] m f 5 [Am) - Am~ f 2 + As + S 2 - 5<r 2 )) 
x T s,a,5,m 2 f,m 2 f ,m 2 ? , 

Any symmetric or antisymmetric coupling function is denned F^jj 
= — Fjj. In the above expressions we have used 



(3.340) 



(3.341) 



Fij + Fji or 



A ab 

%ab . 
A ±i 3 ■ 



B 



ab 



B 



ab 
±ij 



_ s~ia /~ib _i_ s~ia r^tb 

- c a r> h + n a c b 

n« . ,n b .. + n a ..n b 

+13+13 -13-13 > 

- c a n b -+ n a c b 

—ij —ij — ii^— ii > 



13 -V 



r-iab 



+ij +ji ^—ij^ — ji 1 



13 -31- 



ab 

-ij 



D 



ab 
±ij 



f\ab 



- c a n b 4- n a c b 

na . nb ..4. r> a . .r> b .. 

^+■13-^+31 - n -13 -31 ' 



(3.342) 
(3.343) 
(3.344) 
(3.345) 
(3.346) 
(3.347) 
(3.348) 
(3.349) 
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Ch, = A2 1 (A«^)* ± A«« (Ai{f )* , (3.350) 
«i« = Ai'i 1 (Ajfi" ) * ± A!/ a f (A<^) * (3.351) 

4. Conclusions 

We have presented exact cross section expressions for all tree level two body final state 
processes XiXj ~^ A. No approximations have been made in these calculations except for 
ignoring most CP-violating couplings in the SUSY sector. These calculations are useful 
in the proper determination of the neutralino relic density, specifically in the application 
of any constraint on the dark matter density similar to 0.1 < Q x h 2 < 0.3. The dark 
matter density constraint can be used to rule out large portions of parameter space for 
many models of broken super symmetry. Our results encompass the standard calculation, 
in which both initial state particles are the LSP (i = j = 1). Furthermore, our results 
also allow inclusion of all tree level contributions from coannihilation at least one of the 
initial state particles is not the LSP (i and/or j / 1). Such coannihilation contributions 
are important in mSUGRA models when the higgsino component of the lightest neutralino 
is significant and are also important in string-derived models with non-universal gaugino 
masses. 
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